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AÇO TOR 


O AÇO DE ALTA RESISTÊNCIA 
PARA BETÃO ARMADO COM 
MAIOR VENDA NO MUNDO 


TENSÃO DE SEGURANÇA 2 400 Ke /cm 2 
SEM GANCHOS 


Especifique 
AÇO TOR 
com o símbolo 
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G. D. U. 378.245 


A PROPÓSITO DO PRIMEIRO DOUTORAMENTO 
NO INSTITUTO SUPERIOR TECNICO 


NOÉMIO MACIAS MARQUES 


Professor de Física na Escola Naval 


No final do ano lectivo transacto, realizou-se o primeiro doutoramento no Instituto Superior 
Técnico — doutoramento em Engenharia Electrotécnica — e este facto não pode naturalmente deixar 
de ser apontado como dos mais relevantes da nossa vida universitária. Para sublinhar o seu signifi- 
cado ocorreu-nos pedir o acolhimento da Revista Técnica e a esclarecida atenção dos seus leitores. 

A realçar aspectos essenciais do nosso breve comentário, julgâmos de interesse inserir alguns 
excertos do preâmbulo à argumentação do Professor António da Silveira, a quem agradecemos a 
gentileza de nos ter cedido e retocado as respectivas notas, expressamente solicitadas para este efeito. 


£ E E 


A reforma do plano de estudos de Engenharia, estabelecida por diploma de 1955, fundamentou-se 
nas conclusões dós debates havidos em torno de experiências pedagógicas de diversas origens, para 
firmar uma orientação, aliás Já antes seguida por alguns mestres entre nós: o ensino de base deve ser 
constituído por matérias de ciência pura ministradas, com o devido rigor, em extensão e profundi- 
dade. Pressupõe-se que se há-de estruturar o ensino destas matérias sobre um desenvolvimento teó- 
rico de longo alcance e permanente actualidade e que ele deve englobar uma ampla gama de 
elucidativos trabalhos de laboratório. 

Isto pode considerar-se o ponto de partida para conferir em definitivo um carácter eminente- 
mente científico à formação do Engenheiro, quer pela indole de uma parte substancial do curso, quer 
pela sua necessária repercussão em toda a ulterior aquisição de conhecimentos e aprendizagem técnica. 

Como reforço de tal orientação é lícito interpretar ainda a possibilidade, preconizada no mesmo 
diploma, de realização de estágios regulamentares em centros de investigação científica integrados 
na Universidade. Se é certo que, ao nível estritamente escolar do curso, uma sólida base científica 
deve fundamentar a cultura técnica, igualmente certo é que a via mais adequada para atingir a for- 
mação que nestes moldes se esboça, consiste em empreender, ao nível complementar do curso, uma 
iniciação na investigação científica em estágio num Instituto universitário. 

Parece-nos dispensável encarecer o interesse geral destas medidas, numa época em que, por 
motivos bem conhecidos e correspondendo a necessidades bem patenteadas, cada vez mais as activi- 
dades humanas se radicam na investigação científica. Importa-nos aqui sobretudo pôr em relevo 
o papel que uma iniciação científica pode desempenhar na formação do indivíduo. Por sua própria 
natureza a iniciação científica implica a mobilização constante, e a mais exigente, dos recursos 
intelectuais ; oferece horizontes largos e campo vasto à iniciativa criadora ; propórciona oportunidade 
e assunto para ampla e franca discussão das ideias. Por estas vias, em condições adequadas 
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e ambiente propício, fomenta no indivíduo o pleno desenvolvimento da personalidade e o anima 
numa integração social válida. 

Ora é nestas directrizes de reforma — e côntribuindo, a nosso ver, para acentuar o seu real 
sentido — que se enquadra naturalmente a instituição do doutoramento ém Engenharia no Instituto 
Superior Técnico, onde não existia até 1955. 

Na verdade, o trabalho de doutoramento, quer dizer, a elaboração e sustentação de uma tese 
positiva e original — que represente um avanço sensível nalgum domínio do conhecimento — constitui, 
mediante a exaustiva contestação de especialistas autorizados, a prova de maturidade do investigador 
científico. Por isso, o doutoramento dá, em todos os países, em todas as universidades, a garantia 


de uma formação científica autêntica. 
* + * 


O primeiro a apresentar-se a provas de doutoramento no Instituto Superior Técnico foi o 
Engenheiro Manuel Alves Marques, cujo curriculum vitae é um exemplo singular de estrénua dedica- 
ção a uma carreira exclusiva de investigação e ensino. 

Formado no 1. S. T., concluiu o Curso de Engenharia Electrotécnica em 1953, com a classifi- 
cação de 17 valores e foi contemplado com o prémio Bandeira de Melo, Ainda estudante, encetou a 
sua actividade docente no I. S. T., onde tem exercido as funções de Assistente de Física, e foi encar- 
regado durante dois anos da regência da cadeira de Física Geral. 

Bolseiro do Instituto de Alta Cultura desde 1956, estuda, sob a orientação do Professor António 
da Silveira, a estrutura da matéria, com trabalhos experimentais que realiza no Laboratório de Física 
do 1.5. T.. Neste Laboratório encontrou o Engenheiro Alves Marques condições materiais apropriadas 
para o seu estudo por meio de espectrografia do efeito de Raman. Teve no entanto a ocasião e o 
mérito de fazer construir dispositivos complementares da montagem existente, para ampliar e refinar 
a capacidade de observação. Conseguiu assim dar frutuosa continuidade à linha das investigações 
em que se colocou e ao desenvolvimento das ideias para que foi orientado. Sobre os resultados obtidos 
concebeu e elaborou a sua tese de doutoramento, «4 estrutura dos catiões complexos Mg"* (O Ha)s 
e ANt(OHo) em solução», publicada na Revista da Faculdade de Ciências de Lisboa. 

Colaborou em comunicações feitas à Academia das Ciências de Lisboa (Dezembro de 1959) 
e à Academia das Ciências de Paris (Junho de 1961) em cujos Comptes Reniues foram publicados 
em parte os resultados experimentais da sua tese, numa Nota apresentada pelo Prof. Francis Perrin. 
Dá ainda a sua colaboração ao Seminário de Física do Instituto de Alta Cultura em que participou 
com uma conferência proferida em Abril de 1961 sobre «Problemas de estrutura das soluções», 

Em 1960 foi-lhe concedida uma bolsa de estudo para investigar a estrutura da matéria por 
meio dos raios X, mercê de um acordo entre o Instituto de Alta Cultura e a Shell Portuguesa. 


»* * 


A objectividade e a justa medida que devem informar as rubricas de um curriculum vitae 
raramente deixam transparecer alguma coisa da história humana que o enquadra — e essa é por vezes 
suficientemente edificante para ser evocada como lição de vida. 

Do que foi a iniciação do Engenheiro Alves Marques e em que circunstâncias se desenrolou; 
sob que solicitações exteriores venceu, na vontade inabalável em singrar pelo caminho árduo de uma 
autêntica formação científica — disse-o com perfeito conhecimento de causa o Professor António da 


Silveira no preâmbulo da sua argumentação nas provas de doutoramento. Reproduzimos aqui as 
suas palavras textuais: 


«Nestas provas académicas o meu papel deveria ser, naturalmente, mais um papel de rela- 
tor do que de argumentador. Reservo a argumentação para o fim. Primeiro é a título de relutor 
que tenho duas palavras a dizer. 

E natural que os que correram a estas provas — não só estes, mas sobretudo estes — se 
tenham inquirido de como e por que artes o Engº Alves Marques se encontra a fazer um 
doutoramento — o primeiro doutoramento do 1. S.T.. 
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O Eng.º Alves Marques foi um estudante assíduo, procurando cumprir, um estudante uni- 
formemente bem classificado: um «bom estudante» como se diz entre nós, 

Isto, entre nós, são geralmente "simples atributos», não chegam a ser «virtudes», Estes atri- 
bulos, entre nós, geralmente, pouco ou nada abonam em favor do que porventura possa haver 
de autêntico num indivíduo, como virtudes que o imponham acima do comum. 

Não é uma simples impressão, é um juízo formado que poderia ser largamente exemplificado. 

Ainda estudante, entrou para assistente de Física, o que porventura já poderia abonar 
alguma coisa. 

Foi nessa qualidade que um dia me procurou. «Andava eu no 3.º ano» da Faculdade de Ciên- 
cias, quando o Eng. Alves Marques me veio dizer que gostava de se dedicar à Física, de se ini- 
ciar nos trabalhos de investigação científica. 

Devo dizer que logo, logo, não o tomei muito a sério: é que eu tenho uma certa experiência de 
um certo género de «meninos» em busca de qualquer coisa «de fino e de fácil», que querem fazer 
investigação — como dizem — mas que desistem a breve trecho, logo que se apercebem de que inves- 
ligar se escreve com um V... 

Não fiz, bois, grande caso; mas o nosso Engenheiro — que é persistente — voltou e tornou a 
voltar, até que eu resolvi averiguar das suas inclinações *engenheirais». 

E ouvi isto: Que tinha primeiro rejeitado um lugar que lhe haviam oferecido num estabe- 
lecimento do Estado onde fizera um estágio. Depois, que lhe ofereceram um emprego numa com- 
panhia de cimentos. E o proponente, muito ancho da sua condição, acrescentava: não se esqueça 
de que terá casa, água e luz! 

Queriam pôr o nosso Engenheiro a fazer cimento, mas ele — se o não disse pensou-o 
— achou que era cedo de mais para «fazer pó»; e recusou. 

Enfim, mais uma: Um professor do 1.5, 7T. que não está aqui presente perguntou-me um 
dia: Que espécie de indivíauo é que V. tem lá como assistente a quem eu ofereci uma boa situação 
e que recusou alegando desejar primeiro dedicar-se à investigação, adquirir uma formação 
de investigador ? . 

Não sei já o que respondi, mas com estas informações fiquei suspensivo. 

Caspité ! Este engenheiro parece ter um grão de loucura, talvez se passa salvar, 

Tudo isto me fez pensar numa certa obra de Rodin — o imaginário — talvez a mais bela 
e que nunca deixou de estar exposta. Tive a noção de que o Engenheiro Alves Marques sentia 
uma necessidade irreprimível de se libertar, uma ânsia ineluctável de se escapar — de escapar 
à condição humana da generalidade dos colegas do seu tempo. 

Aceitei-o portanto. 

Começou a trabalhar no Centro de Estudos de Física da Faculdade de Ciências de Lisboa. 
Como tinha sido um bom estudante a sua iniciação devia ser acidentada. É que o Eng.º Alves Mar- 
ques tinha de começar por dar provas perante si proprio antes de as dar perante os outros. Como 
tinha sido um bom estudante — um «desalmado» —a sua iniciação devia ser tumultuosa, e foi, 

É não deixou de haver, certamente, quem se empenhasse em duvidar de que ele chegasse um 
dia a fazer alguma coisa. À estes o nosso Engenheiro respondeu sempre, obstinadamente, com 
o silêncio—a um tempo a mais sensível e a mais sublime das respostas, 

Enfim, passados uns dois anos estava a trabalhar com regularidade e muito convencido.» 


“ * 


As provas de doutoramento decorreram perante um júri constituído por todos os professores 
catedráticos do Instituto Superior Técnico e pelos professores catedráticos de Física das Faculdades 
de Ciências de Lisboa, Porto e Coimbra — e com numerosa assistência interessada. 

Tiveram início no dia 16 de Julho com a prova prática, que consistiu na realização de um tra- 
balho de laboratório. No dia 21 de Julho efectuaram-se os interrogatórios sobre os pontos tirados 
à sorte com 48 h. de antecedência, de entre os doze passados pelo Júri e para o estudo dos quais o 
candidato dispõe dos noventa dias que imediatamente precedem esta segunda prova. Decidira o 
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Conselho Escolar do 1. S. T. que fossem os professores de Física presentes no júri a dar os enun- 
ciados destes pontos. Surgiram assim os seguintes assuntos: 


Estatística do radiamento isotérmico. Leis de Rayleigh, de Wien e de Planck. 
Relações de incerteza. Princípio de incerteza. Ilustrações. 
Teoria dos calores especificos. 

Propriedades magnéticas da matéria, 

Condução do calor ; teoria e exemplos da sua aplicação. 
Ondas electromagnéticas nos meios condutores. 

Mecânica dos meios contínuos; deformações e tensões. 

O átomo de hidrogénio em Mecânica Quântica. 

Método de Kron de resolução de redes lineares. 

Reflexão e refracção vitrea e metálica. 

Potenciais electromagnéticos. O campo de uma partícula. 
Fissão e fusão nucleares. 


Do sorteio resultou ser o candidato interrogado pelo Professor Coutinho Braga, da Universidade 
de Porto, sobre «Propriedades magnéticas da matéria» e pelo Professor Almeida Santos, da Uni- 
versidade de Coimbra, sobre «Ondas electromagnéticas nos meios condutores». 

Enquanto que a estas provas se deve atribuir um carácter de exame, entre nós talvez impres- 
cindível, sobre os conhecimentos gerais do candidato ou a capacidade de os aprofundar em curto 
prazo — não resta dúvida de que é a defesa de tese a prova fulcral e decisiva, que polariza todas as 
atenções, na clarividência do exacto conteúdo de um doutoramento. Ela teve lugar no dia 22 de 
Julho, côm argumentações a cargo dos Professores Amaro Monteiro, da Faculdade de Ciências de 
Lisboa e António da Silveira, do Instituto Superior Técnico. 

Através de todas as provas o candidato confirmou as suas excepcionais qualidades e revelou 
conhecimentos profundos, a par de um perfeito domínio dos problemas da sua especialidade. A dis- 
cussão da tese tornou patente a extraordinária importância deste trabalho no campo da estrutura 
das soluções e a valiosa contribuição que fornece para a Física do estado líquido. 

O júri conferiu o doutoramento ao Engenheiro Manuel Alves Marques, aprovando-o com a 
classificação de 7ó valores, muito bom com distinção. 


*o so * 


Uma instituição criada por via legislativa só passa a ter efectiva existência quando os homens 
persistem em dar-lhe realização concreta. A responsabilidade grave do impulso inicial à instituição 
do doutoramento no Instituto Superior Técnico recaiu no Engenheiro Manuel Alves Marques, que se 
mostrou à altura de assumi-la com êxito. Possa o exemplo do Doutor Manuel Alves Marques e a sua 
actividade no Laboratório de Física, onde já foi investido nas funções de 1.º Assistente e prossegue a 
sua carreira científica, tornar-se forte estímulo encorajador dos jovens engenheiros para a reali- 
zação de novos doutoramentos. 

É essa esperança de continuidade e projecção futura que nos impele a evocar aqui, como nota 
final, a fala do Professor António da Silveira, no preâmbulo da sua argumentação, ao salientar a pre- 
sença, no júri, de representantes das três Universidades do País. Dirigindo-se especialmente a esses 
doutos professores, frisou; 


«Lstão aqui os seis físicos da Universidade Clássica mais o físico da Universidade Técnica ; 
são — em linguagem um pouco baroca— os sete oficiais da Física nacional fara talhar, dirão uns 
a capa e balina, dirão outros o «fato de ganga» do primeiro doutor pelo Instituto Superior Técnico. 

São na verdade estas provas de doutoramento as primeiras que se realizam no |. S.T.e 
sendo estas provas académicas as que caracterizam propriamente a Universidade, pode dizer-se 
com propriedade que é pela mão de V. Ex. que o 1. S. T. entra na Universidade «jure et facto». 
Agradeço pois pela minha parte muito sinceramente a V. Ex.“* a colaboração que nos vieram dar 
ajudando a voltar a página e a abrir um novo capítulo na história do 1. S. To. 
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TEORIA DA INFORMAÇÃO 
SISTEMAS DISCRETOS | 


por MANUEL JOSÉ DE ABREU FARO 


Professor de Telecomunicações do 1.5. T. 
Director do Centro de Estudos de Electrônica 
da C. E. E.N. do 1.A.€. 


O objectivo fundamental das Telecomunicações é transmitir correctamente uma mensagem. 

Normalmente o engenheiro de telecomunicações abstrai-se do conteúdo informativo dos sinais e a sua 
preocupação consiste em transmitir sem distorção. 

Ao sinal atribuem-se, então, dimensões no espaço tempo-frequência e exploram-se fundamentalmente os dois 
seguintes métodos de comunicação: «de divisão na frequências e «de divisão no tempos». 

A Teoria da Informação analisa o processo informativo introduzindo conceitos que o esclarecem e funções que 
medem o conteúdo informativo das mensagens. 

Deve-se à Teoria da Informação um conjunto de teoremas fundamentais que permitem averiguar sobre 
o rendimento dos canais que exploramos, nos indicam o sentido da melhoria e asseguram finalmente limites que 
é inútil tentar ultrapassar. 

A Teoria da Informação foi estabelecida em 1948 por Shannon que enunciou e demonstrou os teoremas 
fundamentais. 

Não é possível expor esta teoria nas cadeiras de Telecomunicações I e Il: há que atender a outras 
necessidades mais urgentes 

No entanto, recorremos com frequência aos seus conceitos e assim achamos de interesse organizar o presente 
trabalho que constituirá um elemento de estudo para os que desejarem aprofundar os seus conhecimentos neste domínio. 

Esta é uma razão. Outra, é a necessidade da existência, em português, de livros e escritos para todos aqueles 
que já abandonaram a escola mas continuam com o desejo de se valorizarem actualizando os seus conhecimentos. 

Procuramos dar a este trabalho uma estrutura lógica e acessível que se baseiam na orientação indicada por 
Woodward e Davies e desenvolvida mais tarde por Middleton. Goliman também introduz a Teoria da Informação 
de modo semelhante. 

Procuramos fundamentalmente esclarecer e distinguir as propriedades ligadas à fonte de informação das 
propriedades ligadas ao canal de informação. 

A maior parte das exposições existentes são omissas ou pouco claras sobre este ponto que é delicado e essencial, 

A bibliografia indicada é aquela que de qualquer modo influenciou a presente exposição e que a completa 
da melhor maneira. 


1 — Sistemas discretos: Conceitos fundamentais 
14 — Introdução 


Quando recebemos uma mensagem comunicando que nos coube o ponto número três dum 
conjunto de dez pontos equiprováveis e aceitamos como boa essa informação comparticipamos num 
processo informativo que é de certo modo delicado e embora inconscientemente arbitramos 
alguma coisa. 

De facto: 


* Mantivemos antes de saber o resultado um certo grau de incerteza. 

* Temos a noção de que a incerteza seria menor se o número de pontos possíveis e equipro- 
váveis fosse dois. A probabilidade do acontecimento seria então de !/: enquanto que no 
exemplo considerado era de 4/w. 
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* Temos a noção de que recebemos alguma coisa: uma informação. 


* Acreditando que nos saiu o ponto número três arbitramos que o portador da mensagem não 
alterou o resultado que se verificou no local onde se realizou o sorteio. Assim procedendo 
o grau de incerteza anula-se. 


* Se, no entanto, pomos dúvida na fidelidade do agente portador da mensagem persistirá um 
certo grau de incerteza. À informação não conseguiu anular o nosso grau de incerteza o que 
significa que recebemos uma quantidade de informação menor de que no caso anterior. 


Considerações desta natureza foram o motivo principal da Teoria da Informação. Antes de 
iniciar uma exposição rigorosa tem interesse destacar o significado que nesta teoria se atribui a cada 
um dos elementos dum sistema de informação e os conceitos fundamentais que associa ao processo 
informativo que nesse sistema decorre. 

Tracemos pois o esquema da fig. 1 que representa um «sistema de informação» reduzido aos 
seus elementos essenciais: 


Um destinatário, em B, ligado por um «canal de informação» a uma «fonte de informação» em A 


nte E ne | 
| hd A | Destinatário | 


Canal de Informação 
Fig. 1 


1.1.1 — Fonte de informação: 


Diremos que se constitui uma «fonte de informação» em A se houver expectativa perante 
o que vai acontecer em A 


1.1.2 — Canal de informação: 


É qualquer forma de ligação que se estabelece com o propósito de comunicar a B o que acon- 
tece em A; no sentido de tornar comum ao destinatário localizado em B o conhecimento do resul- 
tado de uma ocorrência em A. 

Um registo feito em B comparado com outro feito em A permitirá estudar experimentalmente 
o comportamento do canal. 

Na inobservância de qualquer acção perturbadora diremos que o canal tem, nas condições em 
que o exploramos e observamos, um comportamento idêntico ao canal ideal sem ruído. 

Um canal ideal sem ruído não introduz alterações, estabelecendo uma ligação certa entre A e B, 

Se pelo contrário observamos uma acção perturbadora diremos que se trata de «um canal 
com ruído». Nestas condições a ligação já não é certa. 

Se os erros provenientes da acção do ruído forem regidos por uma lei de probabilidades 
diremos que o canal com ruído estabelece uma «ligação probabilística» entre A e B. 


1.1.3 — Destinatário: 


É em última análise a entidade a quem se destina a mensagem. Querendo objectivar ó resultado 
da recepção é necessário que o destinatário conheça as características essenciais da «fonte de infor- 
mação» a que está ligado e do «canal de informação» que utiliza. 
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1.1.4 — Grau de incerteza: 


Do exemplo que demos ressalta o conceito de incerteza. 
Introduzir uma função capaz de a medir é um dos objectivos da Teoria da Informação. 
Designaremos essa função por «grau de incerteza». 


1.1.5 — Quantidade de informação : 


No exemplo considerado associamos à informação recebida uma medida que designamos por 
«quantidade de informação». 

Relacionamos ainda a quantidade de informação ao grau de incerteza que persiste após a 
comunicação. 

Resulta desta atitude que na definição de quantidade de informação se fará intervir o «grau 
de incerteza». 


1.1 6 — Lei de probabilidades : 


As fontes físicas de informação e os canais de informação que as leis físicas nos sugerem são 
caracterizados por propriedades estatísticas. Sendo assim, a Teoria das Probabilidades é o instru- 
mento matemático adequado para o desenvolvimento analítico da Teoria da Informação. 

Nesta conformidade admitiremos que: 


a) As mensagens que provêm de uma dada fonte de informação ocorrem segundo uma lei de 
probabilidades que é característica essencial da fonte; 
b) Um canal de informação com ruído pode ser descrito por uma lei de probabilidades. 


Isto significa que se conhece a probabilidade com que uma dada mensagem é transmitida 
correctamente e também a probabilidade associada a determinado erro. 


1.1.7 — Caso discreto: 


Neste capítulo vamos considerar apenas fontes de informação do tipo discreto e finito. 
Isto significa que as mensagens possíveis de uma dada fonte de informação são em número de m 


Mi Me ... Mm 
e que a fonte de informação gera em determinada época uma sequência 
Ms Ma Ms Mi Mi Me Ma Mo Mz.... 
de mensagens cuja ocorrência é regida por uma lei de probabilidades 


1.1.8 — Lei de correspondência : 


Em A fazemos corresponder à mensagem 


Mk 
o simbolo 
Xk 


que transmitido correctamente pelo «canal de informação» seria recebido em B sob a forma 
do simbolo 
Yk 
A acção do ruído poderá ocasionar que em vez de receber y, se receba um símbolo 
Yj 
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Dada a correspondência 
Mk “e Ml 


podemos identificar o símbolo x com a mensagem que simboliza e caracterizar a fonte de informação 
pela sequência de símbolos x; que origina. 


Não havendo ruído a cada 
x: em À 
corresponde com certeza 
y; em B 


e deste modo aos m símbolos transmitidos 


X (xi xs ++. Km) 
correspondem m simbolos recebidos 


Y Oy ya Ym): 


Poderá suceder que o ruído gere símbolos novos, o que significa que transmitido um dado 
X) 
se poderá receber em B um dos quaisquer n simbolos 


Yo Ye Ym Ymgd cce Yn n>m 


Consideremos um exemplo: os símbolos recebidos são sinais telegráficos de um código binário 
de cinco unidades em que só são usadas as combinações correspondentes a três marcas e dois espa- 
ços, portanto, em número de 


simbolos diferentes. 
Teremos neste exemplo m= 10 e n=ãa32. 


1.1.9 — Notações. Relações auxiliares. 


A fim de facilitar a exposição vamos agrupar aqui algumas relações fundamentais da Teoria 


A 


das Probabilidades a que teremos de recorrer. Limitamo-nos à sua interpretação. 
Simbologia : 


p (x;)) — probabilidade de ocorrência de x; em À 

p (yx) — probabilidade de ocorrência de ys em B 

p (x;/y«) — probabilidade de ter ocorrido x; em A na base de se ter recebido y« em B 

p (yx/x;) — probabilidade de se receber y; em B na base de se transmitir x, em À 

p (x; yu) — probabilidade do acontecimento: ocorrência de x; em A seguida da ocorrência 
de yr em B 

p (ye x;) — probabilidade do acontecimento: ocorrência de y; em B antecedida da ocorrência 
de x; em À. 


Relações fundamentais 


1— (x; yk) e (yk xi) traduzem o mesmo acontecimento pelo que p (x; yk) = p (yk Xi) 
2—p (xyx) =p (yrxi) =p (xi) p (yr/xi) =p (vi) p (xi/yx) 
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A relação p(x;yk) =p (xi) p (yk/xi) diz-nos o seguinte: 
A probabilidade do acontecimento «ocorrência de x; em A seguida de ocorrência de yz; em B» 
é dada pelo produto da probabilidade de ocorrência de x; em A pela probabilidade de se receber yk 
em B quando em Á se transmite x;. 
Dir-se-ia cousa análoga de 
P (yr xi) =p (yr) p (xi/yu) 
m 
3—ply)=5 p(xiys) 
I= 


Esta relação é evidente se motarmos que os acontecimentos (x; yz) e (x; y:) ij se 
excluem mútuamente. 


Tn 
4— p(x)= 5 p(yexi) 
k=1 
Tem interpretação análoga ao caso anterior 
5 — De (1) (3) e (4) resulta 


P (yr) = 3 P (x) Pp (yr /x;) 
je 


p (x) = zp (vi) p (x;/yx) 


6 — Notemos ainda que 
m 
Zp(xg)=1 
j=1 


Zzp(Oy)=1 
k= 


m 
Fa p(xi/yr) =1 


 P (yu /) = 1 


E — 


n 
b) 
m mn 

Z E: plyex)=1 

1.2 — Caracterização de fontes de informação e de canais de informação 


1.2.1 — Fonte ergódica de tipo discreto elementar 


Fonte de informação de tipo discreto é qualquer sistema que numa dada época selecciona um 
simbolo x; de um conjunto 
X (x x2 Xm) 


de m símbolos possíveis. A cada símbolo x; corresponde uma mensagem. 
Esta operação repetindo-se N vezes gera uma sequência de N símbolos 


X4 Xb XI XI XA «cc. X3 


cs = e 
N simbolos 


A ocorrência casual e sucessiva dos simbolos x; no tempo, a sua interpretação e expectativa 
que se mantém em relação à natureza dos simbolos seleccionados constitui um «processo informativo». 

No modelo mais simples de fonte de informação não intervêm restrições o que significa que 
as ocorrências de x; e x, numa dada sequência constituem acontecimentos independentes qualquer 
que seja a sua natureza ou a época em que sucedem. 
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Admitindo que a fonte de informação é regida por uma lei de probabilidades, uma sequência 
de N símbolos, N suficientemente grande, não se distingue nas suas propriedades estatísticas doutra 
qualquer sequência de N símbolos gerada pela fonte. 

Quando assim for diremos que a fonte é ergódica e que gera um processo informativo ergódico. 

Quando x; e xx ocorrem com independência, ausência de restrições, diremos que se trata de 
uma fonte elementar. 

Uma fonte que satisfaça a todas estas propriedades será designada por: 

Fonte ergódica de tipo discreto elementar 
a qual é caracterizada por: 


* Um conjunto de m símbolos 
X (x x3....xm) 


* Existência de uma distribuição discreta de probabilidades 


P [p (x4) Pp (xa) .... p (xm)] 


* Ausência de restrições 
p (x; xx) = p (x;) p (xx) 


Numa longa sequência de N símbolos, admitindo que a frequência relativa de ocorrência de x; 
coincide com p (x;), x; ocorrerá um número de vezes que é dado por 


N (x;) = p(x;) N 


1.2.2 — Canal de informação 
Não existindo ruído, o canal de informação estabelece entre A e B a seguinte correspondência 
x; (A) —— y; (B) 


Havendo ruído, um símbolo x; pode dar lugar a qualquer dos simbolos 


YLY2...Ym Ym+ti cc. Yn 
Se a acção do ruído for regida por uma lei de probabilidades existem: 


p (yk/x;) — probabilidade de se receber ys na base de se transmitir x; 
p (x;/ys) — probabilidade de ter ocorrido x; na base de se ter recebido y« 


Qualquer destas distribuições caracteriza completamente o canal, 

Resumidamente teremos : 

* Fonte em À Etr p (mi j=1.,..m 

* Recepção em B [y«,Pp (yu)] k=1...n 

* Canal de informação com ruído descrito por p (y. /x;)) ou p(x;/yk) j=1...m k=1...m 


Das relações de (1.1.9) conclui-se que: 


* Conhecida a distribuição de símbolos y« em Be a distribuição p (x; / yx) é possível determinar 
a lei de distribuição dos símbolos x; da fonte de informação. 
Temos de facto: 
n 
p(xy)= 5 p(yi) plxi/yx) 
| 


k= 


TEONIOA 
82 


* Conhecida a distribuição p(y«/x;) é possível prever a partir da distribuição dos x; a 
distribuição dos y« em B. 
Será: 


p (yo) = s p (x) Pp (yx/x;) 
j=1 


Note-se finalmente a seguinte relação [1,1.1.9] existente entre p(yr/x;)) e p(x;/yk) 


Plyr/x))  ply/ye) 


p(yk) p(x;) 
ld — Grau de incerteza 


Designaremos por 
Us (x;) 
o grau de incerteza que «a priori» se mantém sobre a ocorrência de x; em A, antes de qualquer 
observação em B. 
Designaremos por U (xj/yk) 
o grau de incerteza que «a posteriori» se mantém sobre a ocorrência de x; em A, depois da 
observação de yr em B. 


Associando o grau de incerteza à probabilidade de ocorrência de um acontecimento 
escreveremos 


UV, (x) = u (p(x;)) 
U (x;/ys) == u (p(x;/ys)) 


a função u determina-se a partir das condições que queiramos impor à função «grau de incerteza». 

Começaremos por exigir que: 

O grau de incerteza sobre um acontecimento composto de dois acontecimentos provenientes 
dé uma fonte ergódica elementar deve ser igual à soma dos graus de incerteza que se mantêm em 
relação a cada um dos acontecimentos de per si. 

Sendo assim 

U (x xe) = U (x) + U (x) 
e ainda 
u(plyx))==u(p(x)) +ou (p(x) 


como se trata de uma fonte ergódica elementar 


póxix:) = p(x) ply) 


o que conduz à seguinte equação funcional 


ulxy) = u(x) + u(y) 


cuja solução é da forma: 
u (x) =K log x + kK, 


Determinada a estrutura de u (x) as constantes K e Ko determinam-se a partir do nível a 
que queiramos referenciar o grau de conhecimento e do sistema de unidades que vamos usar. 

Impondo : 

«+ «O grau de incerteza sobre um acontecimento certo deve ser nulo» 


resulta que 
uti) = 6 


uma vez que sendo x; certo p(x;)) = 1 


u(l)=-0O implica Ko =0. 
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* O grau de incerteza sobre um acontecimento de probabilidade 


1 


mn 


representa uma unidade [1], de grau de incerteza. 
Isto significa que devemos ter 


K log — = 1.[1)h 


o que implica 
É 


log n 
Teremos, portanto, 
u(xg)=—-—— Jlogx [1 

Como 

log X = logn x 

log n 

será finalmente 
u(x)=—logax [1]a 
Em unidades binárias teremos 
u (x) == — loga x [1]: == — logs x unidade binária 


A unidade binária designa-se por «binit» (abreviatura de binário digito). 
Usaremos este sistema de unidades. 
» s * 


De acordo com o exposto e com o sistema de unidades que vamos utilizar podemos final- 
mente escrever: 


Us, (x) = — loga p (x) 
U (x;/yx) = — loga p (x; /yx) 
1,4 — Quantidade de informação recebida 


Será medida pela variação do grau de incerteza 
Qu=U, = U 


U, designa o grau de incerteza que a priori, antes de qualquer observação, se mantém sobre 
a ocorrência de um dado acontecimento. 

U designa o grau de incerteza que a posteriori depois de efectuada a comunicação, persiste 
sobre a ocorrência do acontecimento. 

Isto é, recebido um símbolo y« em B a informação que recebemos sobre a ocorrência de x; 
em A é dada por: 


Q (x;/yi) = Uo (x) — U (x;/ys) 


O (x;/ys) = Quantidade de informacão recebida sobre x; na base de termos recebido ys 

U, (x;)) = Grau de incerteza a priori sobre a ocorrência de x; em A. 

U (x1/ys) = Grau de incerteza a posteriori sobre a ocorrência de x; em A na base de ter- 
mos recebido um yr em B. 


Esta definição de «quantidade de informação recebida» é a que se nos afigura de maior sen- 
tido físico. Diz-nos em que medida o destinatário é informado sobre a ocorrência de determinado 
acontecimento a partir de uma observação feita no local de recepção. 
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INDÚSTRIA NACIONAL 
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FABRICA EM REMODELAÇÃO 


MAQUETE DA FARAICA REMDDELADA 


RE 0) DE LAÇÃO 
O CENTRO DE PROJECTOS GUF LAÇÕES 
ESTA EM POSIÇÃO 


DE TOMAR A SEU CARGO O ESTUDO ECONÓMICO DE 
| INSTALAÇÕES EXISTENTES 
PROPÓR AS REMODELAÇÕES NECESSÁRIAS AO AUMENTO DA RENTABILIDADE 
PROJECTAR ESSAS REMODELAÇÕES, QUER FORNECENDO O PROJECTO, QUER ENCARRE- 
GANDO-SE DA AQUISIÇÃO DO MATERIAL E DA EXECUÇÃO DA RESPECTIVA MONTAGEM 


O CONJUNTO DE TÉCNICOS DE OQUE DISPÕE O APOIO DA LONGA EXPERIÊNCIA INDUSTRIAL DA COMPANHIA UNIÃO FABRIL 
NOS VÁRIOS RAMOS DA INDUSTRIA E AINDA À COLABORAÇÃO COM FIRMAS ESTRANGEIRAS GARANTEM A QUALIDADE DAS REALIZAÇÕES 
QUE LHE FOREM ENTREGUES 


O Rato ipi RENO DE PIMDERTOS CA ES 
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FÁBRICA NACIONAL DE FELTROS INDUSTRIAIS 
SOCIEDADE LIMITADA 
APART.: 9 ESCRITÓRIO E FÁBRICA 


ERR E SE OVAR ESTRADA DE S. JOÃO 
TELEG.: FELTROS PORTUGAL 


FELTROS, FILTROS E TECIDOS INDUSTRIAIS PARA TODOS OS FINS: 


EM PEÇAS, PANOS, MANGAS OU SACOS, 
DE QUALQUER FIBRA 

PASTA PARA PAPEL TECELAGENS 
ESA pal ACABAMENTOS 
FIBROCIMENTO TINTURARIAS 
CURTUMES LAVANDARIAS 
FIAÇÕES DE ALGODÃO ESTAMPARIAS 
LOM UA os FSM 0-8 
INDÚSTRIAS QUÍMICAS DESPOEIRAGEM 


IMPERMEABILIZAÇÕES 


TORNA TUDO ESTANQUE 


Impermeabilização é vedação perfeita de: 


Terraços, Empenas, Algerozes, Clarabóias, 
UIBIO Paredes húmidas, Coberturas, Fendas, 
| a | Juntas de dilatação, etc,, etc, 


Fácil aplicação — Elasticidade excepelonal 
Não endurece nem apodrece 


RO ) FATEX À base de borracha + alumínio 


Reforça a acção impermeabilizante do «Rubio», dando às construções um aspecto agra- 
dável — Protege todos os materiais contra a humidade, a ferrugem e a corrosão 


FLUIDO E MASTIC 


RUBIO E ROOFATEX dois produtos que se completam 


Representante Exclusivo : (Olando emantos 


Rua Palmira, 46-A-B-C — LISBOA-! — Telefs. 830163- 8426 39 — End. teleg. «ORFER» 


A existência de ruído impede o completo esclarecimento do destinatário. A existência de pro- 
priedades estatísticas permite definir quantitativamente e num sentido probabilístico o grau de 
incerteza que persiste após a comunicação. 


E K 


Podemos tomar outra posição averiguando qual a quantidade de informação que em A se 
recebe sobre a ocorrência de y; em B na base de na fonte ocorrer x;. 
Antes da selecção da mensagem o grau de incerteza que em À se mantém sobre a ocorrência 
de y;s em B é medido por 
Us (yr) = — loga p (yx) 
Realizada a selecção da mensagem que se supõe ser M; e a que corresponde o símbolo trans- 
mitido x; o grau de incerteza sobre a ocorrência de y; em B é dado por 


U (yx/x;) = — loga p (yx/x;) 


p (yu /x;) = probabilidade de ocorrência de y« em B na base de se transmitir x;. 
De acordo com a definição a quantidade de informação 


O (yx/xj) 
que se recebe sobre yx após a observação de x; será dada por 


O (y:/x) = Us (yo) —U (ye /x;) 


Note-se que 


O (x;/yi) = O (yx/x;) 


De facto: 
O (x/y) =Us (x) —U (9/yx) = [— loga p (x;)] — [— log 2 p (9/yw)] = log” Es 
| 
O (yw/x) =Us (yi) — U (xw/x;) = [— loga p (yx)] — [— log 2 p (yx/x5)] = loge e 
k 


mas, como se notou já, 
pOs/y)  p(yelx) 
p (x;) P (yr) 


o que demonstra a afirmação. 


Uma interpretação interessante destes resultados é a seguinte: 

Consideremos as ocorrências x; em A como resultados de sucessivas experiências. 
Consideremos as correspondentes ocorrências yr em B como resultados de outras experiências 
Entre os acontecimentos x; e y; há uma ligação probabilística descrita por 


P (yr/x;)) ou p(xi;/yk) 


Q (x;/yx) representa a quantidade de informação que em B, a partir da observação do resultado 
yk , Se recebe sobre a hipotética ocorrência de x; em À 

Q (yk/x;) representa a quantidade de informação que em À, a partir da observação do resultado 
X;, se recebe sobre a hipotética ocorrência de y« em B. 

Em qualquer dos casos a quantidade de informação recebida é dada por: 


[anna + ml =| grau de incerteza que — | grau de incerteza que 
recebida do desdela a détemii. | persiste sobre a hipó- 
nada hipótese antes de tese formulada após a 
qualquer observação | observação de um dado 
” x resultado 
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A igualdade 
O (x;/yu) = Q (yr /x;) significa que em 


duas experiências ligadas probabilisticamente a quantidade de informação que se recebe da experiência: 


Hipótese x; 
Resultado yk 


é igual à quantidade de informação que se colhe da experiência 


Hipótese Vk 
Resultado x; 


Não existindo ligação entre A e B,x; e yk são acontecimentos independentes o que significa 


p (x;/yk) = p (x;) (independente de yx) 

p (y:/x;) = p (y«) (independente de x) 
Nesta hipótese : 
P (x; / yk) a logs P (x;) = 


=| 
a p (x;) p (x) 


PO/x) gos Pl — 


3 0 
p (ys) p (yx) 


O (yx /x;) = log: 


o que significa que da observação de um acontecimento independente de outro não se recebe qualquer 
informação sobre este último. 


Pode suceder que 
p (x;/yr) <p (x) 


Nessas condições 


O (y/yi) <<O 
e como 
U (x,/y1:) = Uo (x) — O (x;/yx) 


U (x;/yx) > Uo (x) 


resulta 


Isto significa o seguinte: 

Receber uma quantidade de informação negativa implica que o grau de incerteza final é maior 
que o inicial. 

A experiência realizada (comunicação efectuada) não esclarece mas pelo contrário torna mais 


duvidosa a hipótese formulada. 
» + » 


Esclareçamos a teoria exposta com um exemplo. 
* — Fonte de informação: 
Supomos que só existem dois símbolos xix)2 m=2 
ES 


p (xi) = Ea p (x2) = 


A 
* Lei de correspondência: 
Não havendo ruído xi exs seriam recebidos sob a forma 
vi Yy2 
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* Canal com ruído: 
Não se criam símbolos novos q = 2 
Probabilidade de transmitir correctamente um xi é “mu. 
Probabilidade de transmitir correctamente um xz é “ly. 
Isto significa que ; 


63 po 60 
Pl) = ep(y2/x2) = 64 
Da relação 


ii] 
3 plye/x) = 1 
] 


K=l 
resulta 
[x) =— (Gipa= 
plya/xi) = = ep(yi/x3) = 
da relação 
im 
p(ys)= 3 p(x,) plye/xy) 
fa) 
obtém-se: 


p(y) = p(xo) plyi/x) + p(x:) plyi/xo) = — a 


P(yo) = p(x) plyz/x1) + p(xa) plya/x3) = - Ea 


Estamos agora em condições de determinar a distribuição p (x;/yx) a partir da relação 


púx/yr) = FAIA p(ye/x) 
p(yk) 
o que conduz a 
63 60 
p(xi/y1) o p(x2/y2) E 
4 1 
p(x2/y1) = Rn p(xi/ya) = Fa 


Podemos agora determinar 


Q (x1/y1) = Quantidade de informação que se recebe sobre x1 na base de termos recebido yi 


e 
Q (x:/y1) = Quantidade de informação que se recebe sobre x: na base de termos recebido y1. 
O (xi/y) = UV (x)—U (x1/y1) = [—logz p (x1)] — [— loga p (x /y1)1. 
Q (x2/y1)=U, (x) —U (x2/y1) = [loga p(x2)] — [— loga p (x2/y1)] 
ou seja 
1 63 a 
O (xi/yi) = [- loga T| — |- loga FE | = (1— 0,09) = 0,91 binit 
1 4 é si 
O (x2/y1) = |- loga | e |- loga =| = (1— 4,06) = — 3,06 binit 


No primeiro caso recebeu-se uma quantidade de informação positiva 0,91. Se não houvesse 
ruido, p(xi/yi) = 1 e U(xi/y1) = O o que significa que se anularia o nosso grau de incerteza e se 
receberia mais informação, 1 binit em vez de 0,91. 
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No segundo caso a quantidade de informação recebida é negativa o que significa que o grau 
de incerteza a posteriori é superior ao grau de incerteza a priori. 

Na realidade, recebido um y1 é pouco provável que tenha ocorrido um xa. 

Podemos determinar : 

O (y1/x1) = Quantidade de informação que se recebe sobre y1 na base de ter ocorrido x, em A 

| 1 6: 63 | 
O (yi/x) =Uo(y) — U (yi/x)=| — logz — —— | — | — logs — | =0,94— 0,03 = 0,91 
2 64 | 64 

ou seja, como se notara, 


O (y1/x1) = Q (x1/y1) 


1.9 — Valor médio da quantidade de informação recebida. Entropia. 
A definição 


| Quantidade de informa- | ec | Grau de incerteza | - | Grau de incerteza | 
ção recebida a priori a posteriori 


conduz ao seguinte: 


«A quantidade de informação recebida é positiva se associarmos o resultado observado à 
hipótese mais provável. Se pelo contrário associarmos esse resultado a uma hipótese altamente 
improvável, a quantidade de informação recebida é negativa-» 

No exemplo que se deu admitamos que o canal não tem ruído. 

Nessas condições teríamos: 


p(x/y)=1 p (x:/y3) =1 
p (x2/y1) =0 p(xi/y3) =0 
e então 
O (x2/y1) = Us (x) — U (x2/y1) = loga P (x2/y1) — 
p (x3) 


Isto está certo mas é incómodo pois que condiciona a informação recebida a uma atitude 
arbitrária. 
De facto a quantidade de informação que se receberia da sequência 


Ya y3 YLYI YO ce 


dependeria do símbolo x; que se quisesse associar a cada um dos símbolos y: observados. 

A maneira de resolver a questão será tomar o valor médio da quantidade de informação 
recebida. 

Na sequência considerada, o número de vezes que se observa a transformação 


A —* VE 
é dado, de acordo com a lei dos grandes números, por 
N 6 yr) =p (x yk) N 


em que 
N = número de simbolos da sequência 


p (x; yu) = probabilidade da transformação x; — yk 


O valor médio da quantidade de informação recebida será então 


2 a N (4 yi) O (x/ yr) 


1=0Q (4/yi) = lim * a ia = 2 > p (x; ys) Q (xy/yx) 
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A transformação (x; — yk) representa o acontecimento “(x; yk) —ocorrência de x em À 
seguida de ocorrência de yx em B” idêntico ao acontecimento “ (yk x) — ocorrência de y« em B 
antecedido da ocorrência de x, em A, 

Como se notou: 


piyy)=—p (9) p(yr/x) =p (ye x) =p (ye) P (x /yi) 
utilizando estas relações vamos dar uma outra forma a I. 
Teremos: 
I=2 é p (x yx) O (x ya) = é p (xy yi) [US (4) — U (xy/ yx) 
j 
Calculemos separadamente cada termo. 
E 2 p (x yr) Vo (4) = — E 2 p (x) p (yr/xy) loga p (x) 
= — e p (xy) loga p (x) 2 p (yr / xi) 
dói - p (x) log p (x) 


Note-se que | 
: P (yr/x) = 1 
Façamos: 
— E p (x;) loga p (x;)) = H (X) 


Calculemos o segundo termo mas usando a relação : 


pixy)=Pp (yex)=p (yo) p Og/yk) 


Vem 
gi Pp (xi yo) U (x/yi) = = p (yr) p (x;/ys) loga p (xy/yx) 
=> pn) —2 p (x5/yi) logs p (4/6). | 
=. e Pp (yo) H(X/yvi) =H(X/Y) 


Em resumo: 


Mir médio da quantidade de informação recebida 
sobre X (x;) pelo conhecimento de Y (yu) 


|=em0=1=H09-H (X/Y) 
H (X) = — a p (x;) loga p (xj) 

H (X/Y) = Pp (yo) HX/ yr) 

H (X/yx) a p (x/yx) loga p (x;/ys) 


Analogamente se obtinha: 


| médio da quantidade de informação recebida 
sobre Y (y«) pelo conhecimento de X (x;) 


|- Q 0x/x) =1=H() —H(Y/X) 
H (9) =— 2 p (yk) log: p (yx) 

HER) = p (x) H(Y /x) 

H(Y/x) =— 2 p (yr/x)) loga p (ye /x)) 

1.5.1 — Entropia: 

As funções H introduzidas são normalmente designadas por : 

H (X) = entropia da fonte de informação em A, 
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H (Y) = entropia da fonte de informação em B (saída do canal). 
H (X/Y) = entropia da fonte condicionada ao conhecimento de Y. 
H (X/Y) = entropia da saída do canal condicionada ao conhecimento de X. 


Estas funções gozam das seguintes propriedades : 


* H(X) >0 

* H(X) é máxima para p(xi)=-p(x)) =p (xm)= e , sendo H (X)mazx = log m 
* 0< H(X/Y) &« H(X) 

*+ 0< H(YX) « H(Y) 


Resulta daqui que 
I=H(00-H(XNY) = HOY) - H(/X) > 0 


1.5.2 — Interpretação dos resultados : 
* Entropia da fonte H(X) 
H (X) = — 2 p(xj) log p (x) = 2 p(xj) Uo (x) 
| ) 


representa o valor médio do grau de incerteza que se mantém a priori sobre as ocorrên- 
cias x; em À, 
Não havendo ruído 
H (Y/X) = H(X/Y) = 0 
e então 
J = HtX) = HO) 
pelo que representa também o valor médio da quantidade de informação que se recebe 
usando um canal sem ruído. 
* Entropia condicionada H (X/Y) representa o valor médio do grau de incerteza que após 
a recepção persiste sobre os acontecimentos em A. 
* H(Y) e H(Y/X) têm interpretações semelhantes às anteriores. 
* Valor médio da quantidade de informação recebida I. 


I=H(X) — H(X/Y) = H(Y) — H(Y/X) > 0 


é agora uma quantidade positiva ou nula o que resulta de se ter tomado o valor médio da 
quantidade de informação elementar 


Q (x;/yx) = [— loga p (x)] — [— loga (x;/y%)] 
O segundo termo de Q (x;/yk) contribui para I com 
— P (x//yx) log p (x;/yx) 
Admitindo que p (x;/yk) = 0 
será — plog:p = lim [—- plog: p] = 0 
p— O 
Embora se mantenha significado para uma quantidade de informação negativa, o valor 
médio da informação é positivo ou nulo, o que resulta de os casos mais improváveis 
serem aqueles que menos ocorrem. 
Continua a ser lícito fazer uma hipótese altamente improvável mas essa atitude é pesada 


no cálculo da informação recebida pela probabilidade da ocorrência que é a mesma da 
hipótese que se formulou. 


TEONIOCA 
90 


+ Desde que o ruido seja tal que 
yk é independente de x; 
o que significa em última análise que 
Pp (x/yu)=p (x) 


p (yr/x))=p (yk) 
resulta 


H(XY) = H(X) 
H(Y/X) = H(N) 
e finalmente 
I== 0 


o que significa que deixando de haver ligação probabilistica entre A e B a informação que se recebe 
é nula. 

Este resultado não se aplica só à média mas resulta como é evidente de quantidades de infor- 
mação recebidas nulas por cada símbolo y« recebido. 


Q (x;/yx) = [— log2p (x;)] — [—logrp (x;/y: )]=0 
1,6 — Unidades 
+ O grau de incerteza mede-se em unidades de incerteza [1]n. Adoptando n=-2 teremos 
uma unidade binária que se designa por «binit». 
* À quantidade de informação, diferença de graus de incerteza, mede-se também em «binit». 
* O valor médio da quantidade de informação I é por definição uma grandeza que multipli- 


cada por N nos dá a quantidade de informação recebida de uma sequência de N símbolos. 
Nesta conformidade a unidade de medida de I será 


[1] = binit símbolo! = binit/símbolo 
2 — Sistemas discretos: Capacidade do Canal de Informação 
24 — Introdução 


Admitamos que não há ruído. O canal transmite correctamente os diversos símbolos obser- 


vando-se a correspondência | 
x; (A) — y; (B) 


A quantidade de informação que em inédia se recebe em B é dada por 


I=H(X)-—-H(X/Y) binit /símbolo 
Não havendo ruído 


HAIA )=D 
e então 
I=H()=— a p (x;) logs p (x;) binit / símbolo 
x kk 


Averiguemos agora sobre o comportamento do canal no tempo. Que quantidade de informação 
se recebe na unidade de tempo ou, o que é o mesmo, que quantidade de informação se transmite na 
unidade de tempo? 

Identificaremos a quantidade de informação transmitida com a quantidade de informação rece- 
bida no sentido de que aquilo que efectivamente se transmite é o que efectivamente se recebe. 

A grandeza que mede a quantidade de informação transmitida na unidade de tempo designa-se 
por R 

[R] = binit [unidade de tempo]! 


Para calcular R é necessário atribuir a cada símbolo uma determinada ocupação de tempo. 
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Isso pode-se fazer dum modo simples e que esquemâticamente se representa na figura 2 


A x; xy M km 
T, T, T, Im 
8 T, E T. Tm 
y ), y Ym 


Tempo de Transmissão 


Fig. 2 


Admite-se que a cada símbolo x; corresponde um intervalo de tempo T, contíguo ao intervalo 
Tx correspondente ao simbolo x: que se segue. Supõe-se ainda que após um determinado tempo de 
transmissão se observa em B uma ocupação do tempo idêntica por símbolos y; y« etc. 

A duração média de um simbolo será dada por 


Am T, p (x) 


A quantidade de informação transmitida na unidade de tempo é dada por: 


o — Q(N) 
NT 


O (N) = quantidade de Informação transmitida (ou recebida) numa sequência de N símbolos 
N = número de simbolos da sequência 
T = duração média de um símbolo. 


O (N) calcula-se facilmente, será 
ON =NI=NH(G) 
por definição de valor médio de quantidade de informação. 


Teremos portanto 
“H (X) so Fi P 6) loga p (x) 


1 2 T, p (x) 
| 


R 


Havendo ruído R será dado por 


q HO H(X/Y) 


É 


Surge aqui uma dificuldade devido ao aparecimento de símbolos novos pelo que a questão não 
será abordada com generalidade. No caso do ruído admitiremos que todos os simbolos têm a mesma 
duração inclusive os novos símbolos. 

Sendo assim: 

==... In=T 
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At) ANVIG 5 


— 10º 


f, HERMOFL UID 


| 


+350º 


“ Aquecimento 
industrial 

a alta temperatura 
e sem pressão”. 


e AE RS 


9, AVENUE HOCHE, PARIS 8º 
Téléph.: MAC-MAHON 47-14 


Somente cio 
AANTLHA NOMEL 


aperta [o] Rene 


a, 


FABRICO 
e ACO 
e ACO CADMIADO 


e ACO CADMIADO | 
BICROMATADO 


e BRONZE FOSFOROSO 
e BRONZE BERYLLIUM 

| metricos : 2a 33 mm 
Fera, poros | polegadas : |/l6 a |” 3/8 


- Qualidade e segurança de apêrto 
- Bons contactos eléctricos 


= Aproveitamento de tempo 
na montagem 


* 
Representante : 
GOMES DE CASTRO & IRMÃO, LDA. 
Sede : Filial : 
Rua Guedes de Azevedo,229 Rua dos Remolares, 12-4º« Sala À 
Telefs. 20.658/25.544/36.120 Telefone 366.583 
PORTO LISBOA 


PUB. G. BAUDEL - PARIS 
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ASEA 
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mentos eléctricos industriais 
da Europa Setentrional 


Representada em “Lortugal pela 
ASEA ELÉCTRICA L.” 


SEDE EM LISBOA: | 
RUA DE ARTILHARIA UM, 104, 4.º-DT. 
TELEFONES: P.P.C. 689017/8/9 


FILIAL NO PORTO: 
| RUA DO CAMPO ALEGRE, 144 - TELEFS.: P.P.C.62106/7 
END. TELEG.: "ASEAPORT" 


“IMPERMEÁBILIS" 


O IMPERMEABILIZANTE OFICIALMENTE 
ACONSELHADO EM TODAS AS ARGA- 
MASSAS DE CIMENTO 


Prédios, empenas, caves. — Depósitos para água, vinhos 

e aguardentes. — Facultam-se certificados de ensaio 

dos Laboratório Nacional de Engenharia Civil e da 
Junta Nacional do Vinho 


Distribuidor geral para todo o Império Português: 


ANSELMO DE MATOS 
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+ + * 
A capacidade do canal de informação será, por definição, o maior valor que R pode tomar. 


2.2— Capacidade do canal sem ruido 
2.2.1 — Sistema sem restrições 


Num sistema sem restrições a ocorrência do símbolo x; é sempre possível qualquer que seja j 
e qualquer que seja o simbolo ou símbolos que o antecedam. 

A capacidade do canal será, por definição, 

C == [Valor máximo de R] binit . seg”! 
A determinação de C implica, portanto, a resolução do seguinte problema: 
Dados : 
Já Ta: Im 
determinar a distribuição 
P (p (x1) p (xa) ... p (Xm)) 

que torna máxima a função 


2 p(x;) logs p (x;) 
Res E Rg 
' T,p (x) 


sujeita à restrição 
EF=2p(x)-1=0 


Utilizando o método dos multiplicadores de Lagrange obtemos m equações. 


9 R ) od F =” 
op (xy) dp (x) 


O j=1..m 


Eliminando à é fácil chegar ao seguinte resultado 
p (4) = W-7 


em que W é a maior raiz da equação 


ZW =1 


T, é uma quantidade que se considera adimensional e dá-nos a medida do intervalo de tempo 
relativamente à unidade de tempo que se escolher. 


Isto é o 
T, = T, [unidade de tempo] 


Substituindo em R os p (x;) obtém-se : 
C= logaW binit seg”! 


A interpretação do resultado obtido para p (x;)) é imediato: os símbolos de maior duração 
devem ser os menos prováveis para boa eficiência do sistema, 


Vamos dar exemplos de aplicação : 
1.º Exemplo: Determinar a capacidade de um canal em que Ti = T=Tm=T 
À equação = 
2W == 1 
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assume a forma 


= 
mW =1 
donde resulta: 
log: TA às di m “ 
Ui 


e finalmente: 


E = a logs m binit seg” 


Conclusão: A capacidade do canal representa a máxima quantidade de informação que é 
possível transmitir na unidade de tempo. 

No caso de usarmos simbolos de igual duração, esse máximo atinge-se quando os simbolos 
forem equiprováveis. 


2.º exemplo: Consideremos outro exemplo : 


Dois símbolos xy e xz Ti=20m seg I2=40m seg 
Unidade de tempo 20 m seg 


Tempos T T=1 T:=2 


Teremos 
NU + Wº=1 W = LEVO go log W = 0,696 = 0,7 
do que resulta : 
C=0,7 binit [unidade de tempo] = 0,7 [20 m seg] 
re binit seg”! = 3535 binit seg! 


p 6x1) = (1,62) * = 0,617 
p (x2) = (1,62)? = 0,383 


2.2.2 — Sistema com restrições 


Limitamo-nos a indicar a natureza das restrições que se consideram e a dar o resultado que 


se obtém para a capacidade do canal (Shannon — The Mathematical Theory of Communication). 
Restrições : 


=— fe) 
T; = duração do símbolo x, permissível no estado «a 
e que leva o sistema ao estado [ (*) 


A capacidade do canal é dada por 
C = logs W binit [unidade de tempo]! 


em que W é a maior raíz da equação 


2 W — dg |=0 
; 


à,3 (Símbolo de Kronecker) 


(*) Admite-se que a fonte de informação é susceptível de assumir vários estados. Cada estado distingue-se pela 
natureza dos simbolos que lhe são permitidos. Em resumo, a consideração e a noção de estado constitui um esquema 
que traduz convenientemente um sistema com restrições. 
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Vamos exemplificar : 


Canal alimentado pelos símbolos telegráficos do código Morse que são: 


Ponto ———— T; = 2 unidades Símbolo x« 
Traço —— —» T» = 4 unidades Símbolo xa 
Espaço letra =—» T; = 3 unidades Simbolo x 
Espaço palavra —— T; = 6 unidades Simbolo x« 


Temos dois estados: 


Espoço Palavra 


Espaço Letra 


Es PT 
A o SD 
Traço 


Ponto 


* No estado 1 é permissível a ocorrência de: 


Um ponto que leva do estado 1 ao estado 1 > Ti! — 
Um traço » » » o » IL» » 1 > Tt=4 
Um espaço letra que leva do estado 1 ao estado 2 — Tj? = 3 
Um espaço palavra que leva do estado 1 ao estado 2 —> Til? = 6 
* No estado 2 é permissível a ocorrência de: 
Um ponto que leva do estado 2 ao estado 1 > Ti! =2 


Um traço » x+ » x 2 » » 1 = To! ms À 


Podemos organizar um quadro referente aos diversos símbolos e mudanças de estado que 
promovem : 


11 12 
ng | | | é — 
Tu! Tot | ==: =s ssgiê | = Ta!? | T,t 
2 | 4 | E a E 3 | 6 
mms 
2 4 ms oi - e ss sos 
21 2% 
É fácil escrever agora a equação em W. 
==8 —4 cu —6 
W+W—1 W+W 
== 0 
— 8 —4 e: ( 


W+W 
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o que conduz à equação: 
10 -2 


W+W: LNEW+A WA W=á ss é 
A maior raiz vale 


W = 1,453 
A capacidade do canal valerá: 
C = logs W = 0,54 binit [unidade de tempo]-' 
2.2.3 — Outra definição da capacidade do canal sem ruido: 
Outra maneira de definir a capacidade do canal na ausência de ruído é a seguinte: 


c= tim 108: N (1) 


binit [unidade de tempo]' 
T+ oco T 


em que N (T) designa o número de sequências que é possível organizar no tempo T com os 
simbolos de que dispomos no canal. 

Esta definição conduz aos mesmos resultados da que anteriormente se deu. o 

De facto o número de diferentes sequências que se podem conceber no tempo T tende para 


is T 
N (D)=AW 
A == constante que depende da natureza dos símbolos e restrições existentes. 
W == tem igual dependência e determina-se a partir das equações consideradas em (2.2.1) 


e 223: 
Substituindo N (T) na expressão de C obtém-se: 
C = logs W 
Como W é o mesmo dos casos anteriores a capacidade é também a mesma, 


23 — Capacidade do canal com ruido 


Tem a mesma definição que no caso do canal sem ruido. 
Neste caso admitimos que todos os símbolos têm a mesma duração, pelo que 


R = - [HO)-H (X/Y)] 


A capacidade será então: 
C — [Máximo valor de R] binit seg” 
Os dados do problema são: 
Número de simbolos 
Propriedades do canal quanto ao ruído — p (yr /x)) ou p (xj/ys) 
Será necessário determinar a distribuição 


p (x) p (x3).. p (Xm) 
que torna R máximo 


3 — Sistemas discretos: Teoremas fundamentais 
3.1 — Códigos. Codificação 


Até aqui considerámos o sistema de informação reduzido a três elementos essenciais : 
Fonte de informação 

Canal de informação 

Destinatário 
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tendo-se utilizado a correspondência 


Mk — %j —————+» Yk —+ Mi 
Canal sem 


z Destinatário 
ruído 


Fonte 


Isto é, a cada mensagem My fez-se corresponder um símbolo xx. 

Cada símbolo xy, na ausência de ruído, origina yu em B. 

O destinatário associa a ys a mensagem Mk. 

Deste modo se processa a transmissão de informação. 

Um sistema assim concebido é rudimentar e exigiria tantos símbolos x; quantas as diferentes 
mensagens seleccionadas pela fonte de informação. 

Não é necessário que assim seja e na realidade raramente o é. 

Poderemos transmitir m mensagens diferentes utilizando um número mais reduzido de símbolos. 

Podemos transmitir o que desejarmos utilizando dois símbolos apenas 


X| € Xg% 


Vejamos um modo de o fazer. 
Seja m = 3 
Estabeleçamos as seguintes correspondências : 


Número de ordem Número de ordem Combinações Combinações de Eca a: 
Mensagens na base 10 na base 2 de símbolos x símbolos telegráficos Sinais telmgrívicos 
Mi 1 01 XI X2 Espaço-Marca “Ea 
Ma 2 10 Xx2 X1 Marca-Espaço | | 
Ma 3 11 x2 X2 Marca-Marca | = 


A cada mensagem faz-se corresponder um número de ordem que escrito na base 2 conduz a 
sequências de algarismos 0 e 1. 

A «O» faz-se corresponder xi, a «l» corresponderá xs. 

A x4 corresponde um espaço e a xz uma marca o que conduz em última análise aos sinais 
telegráficos que são típicos de um sistema síncrono de duas unidades. 


A operação que faz corresponder a cada mensagem uma dada combinação de simbolos do 
canal de informação, símbolos telegráficos no caso da telegrafia, designa-se por codificação. 

A operação inversa será a descodificação. 

Código designa o sistema de símbolos utilizados pelo canal e a correspondência que em par 
ticular se estabeleceu entre as mensagens duma dada fonte de informação e possíveis combinações 
dos simbolos distintos que o canal utiliza. 


3,2 — Esquema fundamental de um sistema de informação discreta 


De acordo com o que dissemos, um sistema de informacão discreta traduz-se pelo seguinte 
esquema fundamental (fig. 3): 


Fonte de | | | | “Canal de Centro de o 
informação || codificação informação || descodificação || “estinatário 
TEONTIOA 
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A fonte gera sequências de mensagens. 


No centro de codificação estabelece-se a correspondência : 
My — Combinação x de (x4 x2... xm), 


admitindo que se usam m <. M símbolos x; 
O centro de codificação gera sequências de símbolos 


YE YI cc YmYm+tLcc Yn 


O centro de descodificação opera na base da correspondência 


Xk —* Yk 
a que corresponde em última análise 


Combinação x de (y1 y3 ... ym) — Mk 


Aparecendo uma combinação que não se convencionou (acção do ruido), o centro de desco- 
dificação rejeita-a, podendo pedir à fonte para a retransmitir. Funcionam nesta base alguns siste- 
mas telegráficos com a possibilidade de detecção de erros. 


d.d-— Grandezas características dum sistema de informação 
São as seguintes: 


* Entropia da fonte de informação: 


M 
H (M) = — 5 p (Mi) logs p (Mi) binit simbolo —! 
i=1 
Os símbolos aqui são os simbolos 
Mi Ma ...... Mm 
de M diferentes mensagens. 
Para o centro de codificação, mensagem é qualquer ocorrência a que se faz corresponder uma 


combinação de símbolos 
(x1 Xº +... Xm) 


transmitidos pelo canal de informação. 

Neste sentido, mensagem poderá ser cada uma das letras dum escrito (transmissão letra a letra). 

Poderá ser uma frase dum conjunto de frases convencionais designadas por um número de 
ordem, etc. 

* Código: 

Caracteriza-se pela correspondência que se estabeleceu entre Mk e uma dada combinação dos 
simbolos x; de que dispõe o canal. 

Por exemplo : 

Me —» Sequência xx Xgp -» Xop Px Ou==1 ... m 


ox Cx... Dk —» são índices característicos da mensagem Mk 


* Canal de informação 

Caracteriza-se por: 

Número de símbolos distintos que utiliza xi x2... Xm 

Duração desses símbolos T Tas Im 

Comportamento em relação ao ruído: descrito por p (yr /x;) ou p(xi/yk«). 
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Velocidade de transmissão (de quantidade de informação). 


m 
2 plxi)logp (x) 
= AA) q fé —— (Canal sem ruído) 
li 3 píx) T 
j= 
R = Hit) — H(X/Y) (Canal com ruido) 


Capacidade do canal 
C = [máximo de RJ 


As unidades de medida de R e € são 
[R] = [C] = binit (unidade de tempo)! 


* Quantidade de informação emitida na unidade de tempo (valor médio). 
H (M) representa o valor médio da quantidade de informação que se receberia se não 


houvesse ruído. 
Dividindo H (M) pela duração média de uma mensagem, mais correctamente pela duração 
média das sequências 


Xa XBpo o e Xp 


que correspondem às diversas mensagens, ou seja: 


M E 
Tu= 5 p(M;) [ Tx, + Trat Eta To ] 
j=1 
Ta, , Tê, , To, duração dos símbolos Xay XB Xoys COM Gk, Be Ok = É: sam 
obtém-se uma grandeza: 
H (M) 


M 


G = binit [unidade de tempo]! 


que nos dá a quantidade de informação que é emitida, em média, na unidade de tempo pelo centro 
de codificação e que depende, portanto, das características da fonte e do código que se utiliza. 

* Valor médio da quantidade de informação recebida pelo destinatário. 

Pode-se exprimir na distribuição 


P (Mi) p M3) ... p(Mm) 
e na distribuição 
p (Mi; (9) 
em que W% designa mensagem recebida e M mensagem que foi seleccionada pela fonte. 
Esta distribuição pode-se determinar da distribuição p (x;/y«) característica do canal. 


Será então: 
I=H(M)-H(M/%) 


em que: 


M 
HO(M/WM)= 5 p(M)H(M/M,) 
[=] 


M 
H(M/M)=— 3 p(Mi/M) loga p (Mi /2) 
i=1 
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3,4 — Teoremas fundamentais da Teoria da Informação — Caso discreto 


Estamos em condições de enunciar os teoremas fundamentais da Teoria da Informação. 
Não demonstraremos esses teoremas limitando-nos a dar a sua interpretação A sua demonstração 
encontra-se no trabalho original de Shannon (The Mathematical Theory of Communication). 


I — Teorema — Canal sem ruído (Shannon) 


Dados: 

Fonte de informação de entropia H (M) binit. /siímbolo 

Canal com a capacidade C binit/unidade de tempo. 

É possível codificar a informação gerada pela fonte de tal modo que podemos transmitir com 


uma velocidade de 


F ama ] simbolos/unidade de tempo 


: é positivo e arbitrariamente pequeno. 
Não é possível transmitir a uma velocidade superior a 


= Simbolos /unidade de tempo. 


Notando que 


1 F E = Lá = = = 
—— => número médio de simbolos transmitidos na unidade de tempo 


Tm 


o teorema diz-nos que se pode ter 


Ass 
Tm HM) 
ou seja 
EM cmd 
Tm 


Para H(M) finito «:H (M) é uma quantidade v da ordem de :, vindo 
G=€-) 
o que significa que G deve estar superiormente limitado por € 


IH — Teorema — Canal com ruído 


Dados: 
Fonte de informação de entropia H(M) binit/símbolo e um centro de codificação gerando 
informação ao ritmo médio de G = dona binit /unidade de tempo. 
M 


Um canal de informação de capacidade C binit /unidade de tempo. 
Se G<. € existe processo de codificar a informação gerada pela fonte com um grau de incer- 


teza final 
H (M/9%) 


arbitrariamente pequeno. 
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Se G>C existe processo de codificar de tal modo que H (M/ºi) é menor que 
H(M) — € Tu +: 


: é uma quantidade positiva arbitrariamente pequena. 
Não existe código que conduza a 


H(M/W) <H(M)-—CTm 
+ + + 
3.4.1 — Interpretação 


Canal sem ruído 


Para que um dado canal possa transmitir por qualquer processo de codificação a informação 
gerada por uma dada fonte, é necessário em última análise que o número de sequências que é possi- 
vel organizar nesse canal num dado intervalo de tempo T (T > c<) não seja excedido pelo número 
de sequências diferentes que durante o mesmo tempo possam ser geradas pela fonte. 

O número de sequências que se podem organizar no canal tende para 


TT e O mo 2 


== loga W 
Numa sequência de N mensagens ocorrem 


N (My) = p(My) N mensagens Mi 
N (Ms) = p (M3) N » Ms 
N (Muy) = p (Mu) N » My 


A probabilidade de uma dada sequência particular, ocorrência independente (fonte elementar) 
é dada por 
N p (Mi) N p (Ms) N p (Mm) 
p = p (Mi) GIN] O aaninega p (My) 
ou seja 
= HM) N 


As diversas sequências são equiprováveis e assim o número de sequências distintas com N 
simbolos será dado por 


H (M) N 
E - (M) 
P 


Designando por Tm a duração média de uma mensagem 
T=N Tu 
e então o número de mensagens distintas no tempo T é dado por 


Nm (T) — 2 0 T/Tm 


Devendo ser o P 
N(T) > Nu(T) 
virá 
Ts EM q 
Tm 
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ou seja 


HM) =6CÊC 

Im 

E este o significado do teorema. 
» o» + 

Vamos dar um exemplo: 
Consideremos uma fonte em que: 

p (Mi) = 7/s 

P (Ms) == 1/g 


e um canal em que dispomos de dois símbolos xi e x3 


u>Ti=1 


x — To = 1 
a capacidade do canal obtém-se de: 


WU + WT'=1 C = loga w = log: 2 = 1 binit/unidade de tempo 


Façamos no centro de codificação a correspondência 


My —» xy 
Ma — xa 
Calculemos a entropia da fonte: 
2 7 Po Lg. I é Eai 
H (M)=— 2 p (M;) logs p (Mj) = — E loga z + = loga a = 0,54 binit/símbolo 
j=1 


Calculemos a duração média de cada simbolo (mensagem) 
Tn=pP(M)T+p(M)T=+.1 derem 


Calculemos G, a quantidade de informação emitida na unidade de tempo pelo centro de 
codificação 
HOM) 


Tm 


G = 0,54 binit/unidade de tempo 


o que significa que não aproveitamos completamente as possibilidades do canal que se traduzem 
pela sua capacidade 
C==1 binit/unidade de tempo 
+ + * 
Usemos outro código (código do tipo Shannon-Fano) 
(Mi Mi) = Mi — xa P(M') = p(My P(M) = a 
ER pa é MA P(M') = p (Mi) p(M) = et 


(MM) =Ms— xxx p(M9)=p(M39 p(MD)= ii 


(Ma Ma) = M 4 — xa xa x2 P(M 4) = p(Ms:) p(M3) = = 
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A O EE SP E SS RD AE = 


As mensagens são agora as quatro sequências distintas que se podem organizar com as 
Mi Ms mensagens iniciais. 


Calculemos 
H (M) = — z p (M'i) logz p (M';) = 1,08 binit / símbolo 
Ni 
Calculemos 
Tu= 1 fe 2 + 3 31,36 
Calculemos 
G' = ARO e q 0,79 binit / unidade de tempo 


Tm 1,36 


Melhoramos a velocidade de transmissão de informação de 0,54 para 0,79 binit / unidade 
de tempo. 
Com um código um pouco mais elaborado podíamos aproximar rapidamente o limite 


1 binit / unidade de tempo 
Não poderemos nunca ultrapassá-lo. 
Canal com ruido 


O Teorema de Shannon para o canal com ruído tem um significado profundo. É no nosso 
entender o resultado mais importante da Teoria da Informação e que a justifica como uma teoria 
de grande valor. 

Em última análise diz-nos que a capacidade do canal representa a quantidade de informação 
«que na unidade de tempo efectivamente se pode transmitir. 

Isto significa que se pode escrever simultâneamente 


is | H(M)—H(M/M 
—" "== 
Tu E >0 


-o que conduz a 
H(M/%) = Tm (u—:) = quantidade positiva arbitráriamente pequena. 


A existência de ruído não impossibilita transmitir com uma frequência de erros arbitrária- 
mente pequena. 

A existência de ruído apenas se traduz numa diminuição da capacidade do canal. 

Embora se tenha um H(X/Y) finito existe um código para o qual H(M/W) é arbitrária- 
mente pequeno 


Quando 
Gs E 
será 


H(M/%) 


TEONICA 
103 


está inferiormente limitado por 
— E: Tn=(G—-C) Tum 


É este o significado do teorema fundamental da Teoria da Informação para um canal com ruido 
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THEORY OF INFORMATION 
DISCRET SYSTEMS 


SUMARY 


The fundamental object of the Telecomunication is to transmit correctly a 
message. 

Generaly the telecomunication engineer does not care about the information car- 
ried on the signals and his problem is to transmit without distortion. 

To the signal we must attribute, then, dimensions in the time-frequency space, 
and we explore mainly the two following methods of comunication: «frequency 
division» and «time division», 

The Theory of Information analises the informative process with introduction of 
concepts that clarify it and functions that measure informative contents in the messages, 

The Theory of Information was established in 1948 by Shannon, who enunciated 
and demonstrated the fundamental theorems. 

The present work is an element of study for those who want to get deeper 
knowledge in this held. 

The author gives an easy and logical structure of the Theory of Information 
based on works of Woodward, Davies, Middleton and Goldman. His main aim is to 
ilustrate and distinguish the proprieties of the information source from the proprieties 
of the information channel, 
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PERDAS DE CARGA CONTÍNUAS 


RUGOSIDADE DE TÚNEIS REVESTIDOS E NÃO REVESTIDOS 


1 — Leis de resistência dos escoamentos 
incompressiíveis. 
1.1 — Parâmetros adimensionais. Secções 
circulares. 


Os parâmetros adimensionais que podem figu- 
rar nas leis que regem os escoamentos de fluidos 
incompressiveis em tubagens rectilineas de secção 
constante são os seguintes: 


— Coeficiente de resistência f= Po 
gUÍ 
— Número de Reynolds Re = EUA = E 
e a k 
— Rugosidade relativa A” 


em que é 


Co— Força resistente tangencial média por 
unidade de área do contorno. 

p — Massa específica do fluido. 

U — Velocidade média. 

A — Grandeza geométrica que caracterize a 
forma da secção (diâmetro D no caso 
de secções circulares). 

u — Viscosidade do fluido. 

v — Viscosidade cinemática do fluido. 

k — Rugosidade (hidráulica) absoluta, 


A perda de energia por unidade de percurso ) 
relaciona-se com Z, e, portanto, com f através 
das expressões 


em que é 
Y — peso específico do fluido. 
S — área da secção do escoamento. 


A. QUINTELA 


Assistente do Instituto Superior Técnico 
Da Hidrotécnica Portuguesa 


P — perímetro molhado. 
R — raio hidraulico. 


No caso de condutas de secção circular, por 
ser R =D/4, resulta a fórmula de Darcy-Weisbach 


ID 
= ; 
==, (1) 


1.2 — Regime laminar. 


A lei adimensional que rege os escoamentos 
em regime laminar — no qual a resistência ao movi- 
mento depende sômente dos efeitos viscosos, e a 
rugosidade portanto não intervém — apresenta-se 
como sendo 


f Fa 5 (Re) É 


Neste tipo de escoamento a força de resistência 
ao movimento, para uma tubagem e um fluido 
dados, é proporcional à primeira potência da 
velocidade, o que obriga a lei anterior a apresen- 
tar-se sob a forma 


sendo Cy uma constante dependente da geome- 
tria da secção. 

Para secções circulares a constante Cy toma o 
valor 64 e da conjugação das expressões (1) e (2) 
resulta a lei de Poiseuille (1) 


fo q fc, (3) 


+ [EM 
Os escoamentos na passagem do regime lami- 
nar para o regime turbulento não são estáveis. 


Para números de Reynolds inferiores a 2500 os 
escoamentos têm lugar sempre em regime laminar, 


(1) Pode chegar-se a esta lei por via unicamente teórica, 
a partir da integração das equações de Navier. 
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tendo-se conseguido manter este regime, em con- 
dições excepcionais de laboratório, até números 
de Reynolds da ordem de 70000. Existe assim 
uma zona de passagem do regime laminar para 
o regime turbulento, a qual se designa por zona 
crítica. 

Os escoamentos de água, na sua quási totali- 
dade, não ocorrem em regime laminar. Com 
efeito, para uma temperatura de água de 15ºC, 
as velocidades acima das quais o escoamento é 
turbulento (Re == 2500) são dadas por 


— 0,0028 


U (U em ms”!; Dem m) 


a que correspondem velocidades de 28 cm s”!, 


e 2,8 mm s”!, para tubos de diâmetros respecti- 
vamente 1 cmelm. 


1.3 — Regime turbulento liso. 


Nos esconmentos turbulentos em tubos lisos (tubos 
de rugosidade praticamente nula como os de vi- 
dro, chumbo e cobre) o coeficiente de resistência 
f será somente função do número de Reynolds 


f = q (Re). 


Para tubos de secção circular, esta função 
toma a forma: 


ReV Ff 
2,51 


= 2 log (4) 


Es 
VE 
estabelecida por Karman e Prandtl a partir de 
considerações teóricas e resultados experimentais 
(fórmula de Karman-Prandtl para os tubos lisos). 

As equações empíricas de Blasius e Nikuradse, 
de conhecimento muito difundido, são de utili- 
zação prática mais cômoda que a equação (4) e 
permitem também, dentro dos respectivos limi- 
tes de validade, o cálculo do valor de f para os 
escoamentos turbulentos em tubos lisós. 

O diagrama da figura 1, designado por «harpa 
de Nikuradse», resume os resultados das expe- 
riências de Nikuradse (publicados em 1933) com 
tubagens circulares de raio ro, lisas e com rugo- 
sidade artificial criada por grãos de areia sensi- 
velmente esféricos, de diâmetro constante k, Ve- 
rifica-se que os tubos rugosos até determinados 
valores do número de Reynolds seguem a lei dos 
tubos lisos. O regime dos escoamentos turbu- 
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lentos no interior de tubos lisos ou de tubos 
rugosos enquanto regidos por aquela lei é desig- 
nado por regime dos tubos lisos ou regime turbulento 
liso. 


| rolk = 15.0 


rolk = 126 


2 Eca eee TO Th = 507 
| | R 
| 
| | 
| 
20º E 
404 109 18 Ro 


Fig. 1 — Harpa de Nikuradse. 


1.4 — Regime turbulento rugoso. 


Quando o número de Reynolds cresce acima 
de determinados limites, a resistência ao escoa- 
mento passa a ser determinada praticamente só 
pela turbulência. Atinge-se o chamado regime tur- 
bulento rugoso ou regime puramente turbulento, em que 
o factor de resistência é independente do número 
de Reynolds e, portanto, (ver figura 1) 


 f=ó (k/D). 


Apoiados nos resultados das experiências de 
Nikuradse e com base em considerações teóri- 
cas, Karman e Prandtl estabeleceram a seguinte 
lei para o escoamento em regime turbulento ru- 
goso em tubagens de secção circular 
1 3,7 
—=— = 2 log ——, (5) 

VE k 

Para os escoamentos turbulentos rugosos pode 

utilizar-se também a fórmula empírica de Chézy 


U=CVRI (6) 
ou as fórmulas empíricas monómias 
U=kRº]'. (7) 


Far-se-á referência especial a estas fórmulas no 
número 2, 


1.5 — Regime turbulento de transição. 


No regime de transição entre escoamentos tur- 
bulentos lisos e turbulentos rugosos, o coefi- 
ciente de resistência f dependerá de Re e k/D: 


E tes (o) (Re, k/D). 


Prandtl obteve uma representação da função 
= 6 (Re, k/D) válida para os regimes turbu- 
lentos lisos, de transição e rugosos marcando os 
resultados das experiências de Nikuradse num 
gráfico (figura 2) em que as ordenadas e as 
abcissas são respectivamente 1//f 2 log Dk 
eReVE: Dk. 


lisos. Atinge-se o regime puramente turbulento, 
em que f é independente de Re, quando as aspe- 
perezas penetram fortemente no núcleo turbu- 
lento e a turbulência atinge intensidade tal que 
as acções dissipativas devidas à viscosidade são 
desprezíveis perante as devidas à turbulên- 
cia. 

Colebrook e White vieram estabelecer o elo de 
ligação entre a teoria da turbulência e as suas 
aplicações práticas. Realizando experiências com 
tubagens de betão, tubagens com revestimento 
asfáltico, tubagens de ferro fundido e de outros 
materiais, e de diâmetros desde a fracção da 
polegada até 20 pés, e baseados em coónsidera- 


7 4 10 200 400 1000 
Re VT 
D/k 


Fig. 2 — Função universal da rugosidade. 


É interessante notar que, para as experiências 
de Nikuradse, o regime de transição tem início 
e fim para valores de Re Vf: DK respectiva- 
mente de 8 e 200. Atendendo a que a espessura 
da película laminar à é 


o 116» | 328D 


“O Vith ReV£ 


resulta que o regime de transição tem início para 
9/k ==4 e termina para k 9 = 6. Isto significa que 
os tubos rugosos se comportam como lisos en- 
quanto as asperezas não influem na turbulência 
do escoamento, ou seja, enquanto estão mergu- 
lhadas com certa folga na película laminar, cuja 
espessura é inversamente proporcional ao nú- 
mero de Reynolds. Logo que este atinge deter- 
minados valores, as asperezas, pertubando a pe- 
lícula laminar, influem na turbulência e a lei dos 
tubos rugosos passa a diferir da lei dos tubos 


ções teóricas, propuseram em 1939 uma lei única 
para os escoamentos turbulentos de transição 
em tubagens comerciais 


k 2,91 ) (8) 


VE “37 D já ReVf 


A fórmula anterior, reduzindo-se às duas fór- 
mulas de Karman-Prandtl (4) e (5) respectiva- 
mente para k D - 0 e Re — o, representa os três 
aspectos dos escoamentos turbulentos. 

Na fórmula de Colebrook-White (ou somente 
fórmula de Colebrook) k representa a rugosidade 
absoluta equivalente da tubagem comercial que 
se define como a rugosidade uniforme, a que, 
para o mesmo diâmetro, corresponderia, segundo 
a fórmula de Karman-Prandtl para os tubos 
rugosos, o mesmo coeficiente Ge resistência. 

No diagrama da figura 2, indica-se a curva cor- 
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respondente à fórmula de Colebrook-W hite a qual 
se separa da lei dos tubos lisos para uma abcissa 
menor que a curva representada por Nikuradse, 
e, contrariamente a esta, é sempre descendente. 

Interpreta-se este facto considerando que o 
efeito da rugosidade na turbulência é menos 
brusco nas tubagens comerciais do que nas de 
rugosidade uniforme, em virtude de as asperesas 
poderem ficar umas cobertas pela película laminar, 
outras dela salientes. 

Note-se que os pontos representativos das expe- 
riências de Colebrook e White com tubagens 
comerciais se dispersam em torno da curva cor- 
respondente, afastando-se por vezes sensivelmente 
dela, ao passo que os pontos relativos às expe- 
riências de Nikuradse se aglomeram junto da 
curva respectiva. 

Não obstante este facto e mesmo o de em 
observações diversas se terem verificado, na zona 
dos escoamentos de transição, resultados discor- 
dantes dos de Colebrook e White, continua a 
aceitar-se a fórmula citada, o que não significa 
que se considere estabelecida em definitivo à lei 
daqueles escoamentos. 

A fórmula de Colebrook parece [9] ser aplicá- 
vel a tubagens de «rugosidade isolada», isto é, 
rugosidade formada por protuberâncias separadas, 
de altura e espaçamento aleatórios, provocando 
cada uma um turbilhão independente que dá a 
sua contribuição para a resistência (tipo de rugo- 
sidade da maioria das tubagens comerciais). Não 
é aplicável, pelo contrário, a tubagens de paredes 
caneladas. 

Se para os regimes turbulentos lisos e rugosos 
podem também ser aplicadas respectivamente as 
leis de Blasius ou Nikuradse e a fórmula de Chézy 
e fórmulas monómias, para o regime turbulento 
de transição terá de recorrer-se, e exclusivamente, 
ao emprego da fórmula de Colebrook. 

No número 3 trata-se da aplicação prática da 
fórmula de Colebrook-White e no número 4 dos 
valores de rugosidade absoluta k com especial 
referência aos que se devem prever nos túneis 
com revestimento de betão e nos túneis não 
revestidos. 


1.6 — Secções não circulares 


A fórmula de Poiseuille (3), para o cálculo das 
perdas de carga continuas em regime laminar e 
secções circulares, modificada pela introdução de 
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um factor N, pode tornar-se aplicável a outras 
secções: 


E U 
“a 


J=N (9) 


Poderá ainda utilizar-se para estas secções a fór- 
mula (2) referida a secções circulares (f Re = 64) 
desde que no cálculo de R. se substitua o diâmetro 
pelo quádruplo do raio hidráulico afectado de 
pU4R ) 


v 


um factor p correctivo (R. == 


Indicam-se a seguir os valores de N, obtidos 
por Boussinesqg [9], bem como os valores de p 
que lhes correspondem: 

— Secção circular (totalidade da secção) :N=2; 

p=l. 

— Secções rectangulares: 


Altura largura 


Em pressão | EM superfície 
1,0 0,50 1,78 1,12 
0,5 0,25 1,94 1,03 
0,333 0,166 2,15 0,93 
0,25 0,125 al 0,88 
0,2 0,1 2,38 0,84 
ó; 0,05 2,65 0,75 


O cálculo de escoamentos turbulentos em con- 
dutas e canais de secção não circular, quando se 
utiliza a fórmula de Colebrook, tem sido reali- 
zado substituindo simplesmente o diâmetro D 
pelo quádruplo do raio hidráulico D = 4R. 

Em 1938, Keulegan obteve expressões teóri- 
cas, complexas, para o cálculo dos escoamentos 
turbulento liso e turbulento rugoso em canais 
rectangulares de largura infinita e canais trape- 
zoidais [9], propondo, depois de algumas simpli- 
ficações, a fórmula para escoamentos turbulentos 
rugosos em canais: 

Rr Pro e 
V£ k 


Os valores de f obtidos por esta expressão e 
os obtidos pela correspondente fórmula de Kar- 


man-Prandtl (5) — em que se faz D = 4R — di- 
ferem entre si de 5 e 10º, quando 2R k toma va- 
lores entre 10 e 10º. 


Atendendo a estas pequenas divergências, tem 
sido prática corrente, para o cálculo dos escoa- 
mentos turbulentos em pressão em secções não 
circulares ou em superfície livre, a utilização — 
com a obtenção de resultados aproximados — da 
fórmula de Colebrook em que se substitui o diâme- 
tro pelo quádruplo do raio hidráulico. Este pro- 
cedimento equivale a admitir que o efeito da 
secção é completamente representado pelo raio 
hidráulico —- como na aplicação das fórmulas 
empíricas — e que o efeito da superfície livre é 
despresável. 

A primeira hipótese parece pouco aceitável se 
se considerar que as fórmulas de resistência de 
Karman-Prandtl resultaram da integração, para 
a totalidade da secção circular, das equações es- 
tabelecidas pelos autores para a distribuição de 
velocidades em regime turbulento naquela sec- 
ção, e que da integração para secções diferentes 
“deverão resultar outras constantes de integração. 

A validade da segunda hipótese implicaria 
a verificação de um comportamento igual, sob o 
ponto de vista hidráulico, dos escoamentos em 
pressão e em superfície livre, desde que a secção 
do primeiro fosse obtida da do segundo por 
acréscimo da secção simétrica em relação à su- 
perfície livre. Tem-se porém observado que nos 
escoamentos em superfície livre a velocidade 
máxima não ocorre à superfície — como corres- 
ponderia à igualdade de comportamento referida 
— mas abaixo desta. 

Muito recentemente, Enrico Marchi [11] indi- 
cou factores correctivos do raio hidráulico a 
introduzir na fórmula de Colebrook, de modo a 
serem considerados os efeitos da forma da sec- 
ção e da superfície livre. 

Por considerações teóricas e apoiado em resul- 
tados de experiências realizadas pelo próprio e 
por outros autores, propôs fórmulas para a dis- 
tribuição de velocidades nos escoamentos turbu- 
lentos lisos e rugosos no interior de condutas de 
diferentes secções. A partir dessas fórmulas 
obteve as seguintes leis de resistência : 


— Regime turbulento liso 


; ar 
dE = A, + Blog Re Vf 


— Regime turbulento rugoso 


em que 
— U(4R) | 


Y 


Re 


O valor do coeficiente B é independente da 
forma da secção e igual a 2,00. Os coeficientes 
A, e Ar dependem da forma da secção e da dis- 
tribuição no contorno da resistência tangencial 
unitária. 

E. Marchi reduziu as expressões anteriores às 
seguintes :; 


— Regime turbulento liso 


-——— =2] 2,51 (10) 
— Regime turbulento rugoso 
1 3,7p (4R) 
—— = 2 log —————""— 11 
VE ni * a 


pelas quais se verifica a possibilidade de traduzir 
o efeito da forma da secção por um factor p 
correctivo do raio hidráulico, próprio do re- 
gime turbulento e independente da natureza das 
paredes. 

Gentilini, anteriormente a E. Marchi, havia 
proposto para o cálculo dos escoamentos turbu- 
lentos pelas fórmulas de Karman-Prandtl, a 
adopção dos valores do factor p obtidos para 
o regime laminar. 

Pelas mesmas considerações teóricas que con- 
duziram à fórmula de Colebrook, obtém-se 


1 | et k 2,51 | 

WE pL3I7(4R) ReVtf 

1 | k 2,51 
—— = — 2 log | - e qd UR) 
VE | a + | 


Esta expressão é idêntica a (8), desde que no 
cálculo do número de Reynolds e da rugosidade 
relativa se substitua o diâmetro pelo quádruplo 
do raio hidráulico afectado do factor p 


(D=4R, e Rep ms (6 Ro) 


com R=pR). 


v 


Para os esconmentos em pressão, E. Marchi, utili- 
zando os resultados indicados por vários autores 
para escoamentos em tubos lisos, obteve os se- 
guintes valores de p: 


TTYONIOA 
109 


— Secção quadrada — 1,00 à 1,17 
— Secção triangular equilátera — 1,30 
— Secção rectangular larga  — 0,84. 

(ou anular estreita) 


Note-se que, em conformidade com as consi- 
derações teóricas de E. Marchi, são diferentes os 
valores do factor p para os regimes laminar ou 
turbulento liso. A aplicabilidade ao regime ru- 
goso dos valores do factor p determinados para 
o regime turbulento liso — que resulta daquelas 
considerações — não foi ainda comprovada expe- 
rimentalmente. 

Para os escoamentos em superfície livre, a partir de 
resultados próprios e dos indicados por outros 
autores para escoamentos em tubos lisos e rugo- 
sos (de rugosidades artificial e comercial), E. Mar- 
chi determinou os valores a seguir indicados para 
o factor p, que neste caso traduz a influência si- 
multânea da forma da secção e da superfície li- 
vre : 


Secção triangular equilátera = 1,25-— 1,30 
Secção triangular recta p=1,15—1,20 
Secção triangular larga p=-1,00 
Secção semi-circular p==0,90 


Secção trapézia (= 1/2 hexágono) p==0,90 — 1,00 
Secção trapézia larga p==0,80 
Secção rectangular de economia p==0,95 


(bo = 2 yo) 

Secção rectangular estreita p=0,90 
(bo = 2/3 Yo) 

Secção rectangular larga p=0,80. 
(bo = 6 yo) 


Nas designações relativas às secções rectangu- 
lares, yo e bo significam respectivamente a altura 
de água e a largura do canal. 

E. Marchi indica que o erro relativo dos valo- 
res que determinou para o factor correctivo é 
inferior a 5 */h. 

Para finalizar, cabe indicar a influência intro- 
duzida pelo factor correctivo no valor do coefi- 
ciente de resistência. O factor correctivo difere 
da unidade no máximo de 0,2, com excepção do 
caso das secções triangulares, pouco frequentes 
na prática. A introdução deste factor no cálculo 
do número de Reynolds ou da rugosidade rela- 
tiva provocará nestas grandezas uma variação 
máxima de 20 */o. 

Nos escoamentos turbulentos lisos, esta variação 
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do número de Reynolds implica uma variação 
do coeficiente de resistência (e da perda de carga 
unitária) de cerca de 6 “/,e 4º/, respectivamente 
para Re = 10* e Re = 10º. 

Nos escoamentos turbulentos rugosos, a mesma 
variação da rugosidade relativa determina varia- 
ções do coeficiente de resistência de cerca de 9º, 
71/y e 4º/yrespectivamente para rugosidades rela- 
tivas de 0,05, 0,01 e 0,0001. 


2 — Fórmulas empíricas para escoamentos 
turbulentos rugosos 


As fórmulas empíricas de uso mais generali- 
zado para o cálculo de escoamentos turbulentos 
rugosos são as de Chézy e de Manning-Strickler. 

A fórmula de Chézy (6) foi estabelecida em 
1775 para o regime uniforme em canais a partir 
da consideração de que a resistência tangencial 
média por unidade de área do contorno seria 
proporcional à energia cinética 


Co=À —— (13) 


e da condição de equilíbrio entre a força resis- 
tente num troço ÀL e a componente do peso na 
direcção paralela ao fundo. 


SAL] = PAL. (14) 


De (13) e (14) resulta 


sa TE Tc yu Ee 
pRJ=C 5 
e 
U=CVEL (9 


A experiência demonstrou que o coeficiente C 
depende, além da natureza do material da tuba- 
gem, da geometria da secção. 

Para o cálculo do coeficiente C dispõe-se, entre 
outras, das fórmulas de 


/ VR 


Basta € = —— (em unidades métricas) 
= VR 
— 100VR | 
Kutter C=-————— (em unidades métricas). 
m + VR 


A fórmula de Manning-Strickler é uma fór- 
mula monómia 


V=kKR* y!! (15) 


e pode resultar da fórmula de Chézy, se for 


C=KR”, 


O coeficiente de velocidade K tem as dimen- 
sões L'!3 T-, 

Os parâmetros y, me K são dependentes da 
natureza do material da tubagem e existe ampla 
informação dos seus valores, [5], [12]. 


UD 


dos para qualquer fluido desde que o movimento 
seja turbulento rugoso. 

Referiu-se em 1.5 que o movimento turbulento 
rugoso tem início para Re VE: D/k = 200, ou, 
como o valor de f no escoamento turbulento 
rugoso é dado por (5), para 


a ;D (og 27 D 
= 2002, E 


Da 


Na figura 3 apresenta-se a função 


3,7 D 
16 
E ) (16) 


UD = 23>=<10"* (2 log 


Ds 
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Em 


107º + 
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Fig. 3 — Início do movimento turbulento rugoso. Escoamento de água a 15º €. 


Estas fórmulas — que se aplicam a secções di- 
ferentes da circular e aos escoamentos em super- 
fície livre— são de emprego prático, em especial 
a fórmula de Manning-Strickler, que se pode 
traduzir por um ábaco muito simples. 

Os valores obtidos por estas fórmulas são váli- 


e no quadro 1 os valores das velocidades para 
diferentes valores do diâmetro e rugosidade abso- 
luta, correspondendo ao início do movimento 
turbulento rugoso para o escoamento de água 
a 15ºC (v == 1,15 x 10-*m?s—!). 

Também se representa na figura 3 e se indica 
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no quadro 1 o limite do regime turbulento rugoso 
indicado por Ackers para a água a 15º C (Uk = 
— 0,928; k em mm e U em ms”), a que se fará 
referência em 3.2. A coincidência dos dois limi- 
tes verifica-se para valores da rugosidade relativa 
muito elevados, da ordem de 0,04. 

Sendo as velocidades dos escoamentos nos 
túneis revestidos dos aproveitamentos hidroeléc- 
tricos, normalmente inferiores a 4,00 ms”! e a 
rugosidade absoluta por vezes inferior a 0,5 mm, 
aqueles escoamentos são, em muitos casos, do 
tipo turbulento de transição, como se conclui 
do exame do quadro 1. 

Note-se que as fórmulas empíricas não apre- 
sentam compatibilidade absoluta com a lei de 
Karman-Prandtl para tubos rugosos, [1]. 


tubos rugosos obrigaria K a depender — além da 
rugosidade k característica do material, como se 
indicou — do diâmetro da conduta. 

No diagrama da figura 4 representa-se a curva 
correspondente a (18); neste diagrama uma recta 


corresponde à relação 


D a (KDU 
K 


sendo 6 o valor de b para o qual K se torna 
independente do diâmetro. 

A recta traçada neste diagrama corresponde, 
[4] e [9], à relação 


k = [as] (k em m) (19) 


também traduzida no diagrama da figura 5. 


QUADRO 1 


Limites inferiores da velocidade U ms”! em escoamentos turbulentos 
rugosos. Água a 15º C. 


| | 
Diâmetro | 


CT —————e— | mu | ——e e o 


Limite Uk = 0,928 


(U ms”!; k mm) 


X Rugosidade absoluta (mm) 

O | (8) Pele ]ale) 

m 

be a— 0,1 0,5 1,0 5 10 100 
L & a nm 
m o | | 
o 0,10 16,5 2,60 1,18 0,17 0,07 — 
a: GR 0,50 20,0 3,30 | 1,44 | 0,24 | 0,11 — 
| = 1,00 21,0 360 | 1,64 | 0,26 | 0,12 | 0,007 
= X 5,00 24,0 4,20 2,00 0,33 0,15 | 0,010 
” 10,00 25,6 4,40 2,10 0,35 0,16 | 0,012 
E 15,00 26,2 4,67 2,20 0,37 0,17 0,013 
pl 


Para a fórmula de Manning-Strickler conside- 
rando (1), resulta, em tubagens da secção cir- 
cular 
884 


= FRDR' 


(17) 


Igualando os valores de f de (17) e de (5), 
resulta 


V127g (2 log 3/7 D/K) 


paes Dis 


(18) 
Verifica-se assim que a compatibilidade das 
fórmulas de Manning e Karman-Prandtl para 
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Os pontos de cruzamento da curva (18) e recta 
(19) correspondem aos valores da rugosidade 
relativa para os quais são compatíveis as fór- 
mulas (5) e (15), que, para outros valores de 
rugosidade relativa, fornecem valores de f dife- 
rentes. 


3 — Aplicação prática da fórmula de Cole- 
brook-White. 


3.1 — Ábacos de Moody e Rouse. 


A fórmula de Colebrook-White apresenta por 
um lado a vantagem da sua aplicabilidade a todo 
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Fig. 4 — Compatibilidade das fórmulas de Karman-Prandtl (regime turbulento 
rugoso) e Manning-Strickler. 


ES 
NESSE E E E O PR DO E 


“a 8 5 k(m) 


Fig. 5 — Relação entre a rugosidade absoluta k e o coeficiente K de Manning-Strickler- 


o campo dos escoamentos turbulentos, por outro da fórmula de Colebrook dos seguintes problemas 
lado, a inconveniência de uma aplicação prática práticos: 
complexa, o que constituiu obstáculo à rápida 


difusão do seu emprego. 1º — Dados U (ou Q), ] e k, determinar D; 
Pode pretender-se a resolução por intermédio 2º — Dados U (ou Q), D ek, determinar J ; 
TEONICA 
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3º — Dados D, ] e k, determinar U (ou O); 
4º — Dados D,U (ou O) e J, determinar k. 


A fórmula de Colebrook escrita sob a forma 


U=-2V2gDJ]log (— 2,91) V 


220) 57555) (20) 


permite a resolução analítica directa do 3º e 4º 
problemas, sendo necessário resolver por tenta- 


constante e a recta É Re = 64 correspondente à 
equação de Poiseuille. Sendo aplicável aos movi- 
mentos laminar e turbulentos, é designado por 
ábaco universal de resistência. 

Permite a resolução directa do 2º e 4º pro- 
blemas e obriga tentativas para a resolução 
do 1º e 3º, 

O ábaco de Rouse (figura 7) é também um ábaco 
universal de resistência; o eixo das ordenadas é 
graduado em 1 WE; os dois eixos das abcissas são 
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Fig. 6 — Ábaco de Moody. 


tivas os restantes. Na zona dos esçcoamentos 
puramente turbulentos, a equação (5), à qual se 
reduz a fórmula de Colebrook, já permite a reso- 
lução analítica directa dos quatro problemas. 

Aquelas resoluções analíticas, mesmo directas, 
obrigam à execução de cálculos laboriosos, tendo 
sido por isso estabelecidos vários processos, quer 
analíticos quer gráficos, para tornar cómoda a 
aplicação da fórmula de Colebrook. 

São citados nesta sequência os ábacos de Moody 
e Rouse, muito divulgados. A bibliografia mais 
recente — [6], [9], [10], [13], [14], [15], 
e [16] — refere progressos na matéria, dos quais 
se consideram como mais felizes os ábacos de 
Ackers e as tabelas de Schulz, de que se tratará 
nos números 3.2 e 3,3, 

No ábaco de Moody (figura 6), cujas abcissas e 
ordénadas são respectivamente os logaritmos 
de Re e f, figuram as curvas f (Re) para k/D 
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Fig. 7 — Ábaco de Rouse. 


graduados em logaritmos de Re e de Re V £ 
— V2g]D'/?2. 
Com o ábaco de Rouse consegue-se a resolução 
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Fig. 8 — Ábaco de Ackers, 


